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A COMPUTER PROGRAM FOR ELECTRON PROBE MICROANALYSIS 
IN THE FIELDS OF METALLURGY AND GEOLOGY 


J. I. Goldstein* 

Metallurgy and Materials Science Department 
Lehigh University 
Bethlehem, Pa. 

and 

P. A. Cornelia 

Laboratory for Theoretical Studies 
Goddard Space Flight Center 
Gre.enbelt, Md. 


ABSTRACT 

A comprehensive computer program has been written which should prove 
useful to investigators in both the fields of metallurgy and geology. The program 
is flexible in both input and output and is relatively easy for an inexperienced 
person to use. It is written in the Fortran language and can be used in remote 
terminal operation which is now available in many laboratories. The program 
is built up of a large number of subroutines which can be easily changed as the 
user sees fit or as new correction schemes become available. At present, the 
following corrections are incorporated; the absorption correction of Philibert- 
Duncumb 1 1 2 recently modified by Heinrich, 3 the atomic number correction of 
Duncumb and Reed, 4 and the fluorescence correction of Reed. 5 K a x-ray radia- 
tion from elements of atomic numbers 3-39 and L a radiation from elements of 
atomic numbers 20-92 are considered. The program operates in either of two 
modes; conversion of raw intensity data to composition using an iterative pro- 
cedure, or conversion of compositional data to expected x-ray intensities. For 
problems in geology involving oxides, calculations are made with oxygen con- 
sidered as a matrix element. Specimens and complex standards each containing 
up to nine elements can be treated and several standards can be used for each * 

element measured in order to check the consistency of the answer. The program 
determines internally whether a fluorescence correction is necessary and also 
calculates the initial intensity ratios. The program does not correct for non- 
normal incidence or continuum excitation. An exhaustive print out of the various 
correction factors can be obtained if desired. 


*Formerly at Laboratory for Theoretical Studies, Goddard Space Flight Center, Greenbelt, Md. 
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The program is rather unique with respect to the method of data handling 
and ease of operation. All the critical data for a given element, such as x-ray 
wave lengths, fluorescence yields, atomic weights, absorption edge jump ratios 
and the necessary factors for the calculation of mass absorption coefficients from 
Heinrich 6 are stored on cards or are placed on file at the computer. These data 
are then used as needed for the elements under consideration in each problem. 
The compositional data for each of the element standards are also stored on 
cardu or at the computer. The input data which specifies the problem of interest 
then requires a minimum of information and can be supplied in anon-rigid format. 



A COMPUTER PROGRAM FOR ELECTRON MICROPROBE ANALYSIS 
IN THE FIELDS OF METALLURGY AND GEOLOGY 


I. INTRODUCTION 

The electron microprobe has been used to solve problems in many fields of 
scientific research such as metallurgy, geology and biology. In the area of 
quantitative microchemistry the results of microprobe analysis have not fullv 
met the expectations of the early pioneers. Schemes which correct the x-ray 
data for the effect of atomic number, absorption and fluorescence have been used 
extensively. However at this time precision errors (Ziebold) 7 and errors in the 
input parameters to the various correction schemes (Heinrich and Yakowitz) 8 
limit the precision and accuracy of the analysis to no better than about 2% relative. 
Comparisons between measured intensity data from well characterized standards 
and calculated intensity ratios (Ziebold and Ogilvie, 9 Goldstein, et al. 10 and 
Beaman 11 ) support this conclusion. 

To avoid using any corrections for microprobe analysis, standard samples 
with similar compositions to the analyzed sample must be obtained. In the vast 
majority of cases this is not practical, and the correction schemes must be used. 
As the number of unknown elements, standard samples and data points desired 
increases, the amount of data handling and calculation becomes excessive. The 
use of computer techniques allows rapid handling of data, accurate calculation 
of correction procedures and assimilation of output in a useful form. With the 
development of remote terminal operation and rapid microprobe readout devices, 
the need for computer programs for quantitative electron microprobe analysis 
has developed. 

About 20 computer programs for electron probe microanalysis are now avail- 
able. Descriptions of the various programs are presently being compiled by 
Beaman 12 and need not be discussed here. In general many of these programs 
are either inflexible, are written for specific computer systems or are only able 
to handle a small number of elements or standards. Since problems originating 
in various disciplines are often run in the same microprobe facility, a program 
was designed that could be used to treat both geologic and metallurgical problems. 
We desired a flexible program built on a modular design with minimal input and 
data handling, compatible with remote terminal operation, and incorporating the 
best correction procedures currently available. 
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This paper describes the program we have developed. The program design, 
calculation schemes and input-output are described in sufficient detail so that 
the computer program can be adopted for use in most laboratories. Several 
sample problems are discussed in order to illustrate the use of the program, 


n. CORRECTION PROCEDURES FOR ELECTRON MICROPROBE 
ANALYSIS 

A. Corrections to Raw Intensity Data 


To obtain an analysis for each element in a sample of interest the following 
measurements are made: (1) the x-ray intensity from the element standard 
taken both before and after the analysis, (2) the x-ray intensity from the back- 
ground standard or the continuum radiation above and below the x-ray peak, and 
(3) the x-ray intensity for each element from the sample of interest. Three 
changes are made in the raw data, so that it is in a suitable form for calcula- 
tions of the x-ray corrections. 

1. Detector Dead Time - Detector dead-time causes a discrepancy between 
the true count rate of x-rays incident on a detector and the measured count rate. 
This relation is expressed as: 


N s n/(l - r» r) (1) 

where N is the true count rate, n is the measured count rate and r is the de- 
tector dead time. This correction is applied separately for each element. 

2. Drift - Drift is the slow change in x-ray signal caused by instabilities in 
the electronic system of the microanalyzer. Normally this change is less than 
1% unless unusually long analysis times are used. In this program, if the stand- 
ard intensity for each element differs, from the beginning to the end of an anal- 
ysis, by more than 6 times the square root of the average standard intensity, a 
drift correction is applied. It has been found in practice that if the change in a 
monitored signal is no more than 5% per hour, then it may be assumed that the 
drift is linear in time. A linear drift correction is therefore applied when 
necessary. 

3. Backgro und - Background intensity arises from the continuous x-ray 
spectrum generated in the sample. It is obtained for each element by either 
measurement of the element intensity from a background standard or at spec- 
trometer settings away from the Bragg angle. When samples of interest contain 
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many elements it is difficult to obtain an ideal background standard, that is, one 
which contains all the elements except the element being measured. Therefore 
the background is oommonly measured by the second method. In the computer 
program, the background intensity is subtracted directly after the data is cor- 
rected for dead time. 


B. Correction Procedures for the Conversion of X-Ray Intensity 
to Composition 

It was first demonstrated by Castaing^ that the relationship between the 
mass concentration, C A , of an element, A, in a specimen analyzed in the micro- 
probe and the ratio (k A s I A /I (A) ) of the intensity of the characteristic x-ray 
line of element A measured from the specimen, I A , to that measured from a 
standard, 1(A), is, to a first approximation, linear: 

k A s C a /C(A) (2) 

where C(A) is the amount of A present in the standard. If the standard consists 
of pure A, then 

k AS C A . (3) 

Departures from the simple law of proportionality expressed in Equations 2, 3 
can be large. Therefore if quantitative microchemical analysis is to be attempted 
with the microprobe, it is essential to apply a correction factor: 

C. (correction factor^ . 

k A = I a /I(A) = -Ai- — f - (4) 

C(A) (correction factor) std 

This factor, is needed for the standard when it is not a pure element. The cor- 
rection factor is placed on the right hand side of the equation because it is itself 
a function of C A , C(A) and of the concentrations of the other elements present. 

Following Castaing, most authors have divided the total correction into 3 
separate factors: atomic number, absorption and fluorescence. The atomic number 
correction takes into account the effect of composition or atomic number on the 
intensity of characteristic x-rays generated within a specimen. The absorption 
correction allows for the absorption of characteristic x-rays in emerging from 
the specimen. The fluorescence correction allows for the additional intensity 
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contributed by fluorescence radiation which is due to excitation by x-rays from 
other elements in the specimen. The measured intensity ratio is therefore given 
by; 


X A C A (At.#L [Abe] [flu or] 

k £- — » — (5) 

A 1(A) C(A) [At.#],, [Abs] B [fluorlg 

where u represents the specimen of interest and s represents the standard. In 
most geologic problems C(A) is not equal to 1.0 and Equation 5 must be applied. 

In the following sections the three corrections which are applied to unknown 
and standard Intensities are described in detail. They were chosen as the best 
correction schemes available to date. It should be noted that none of these 
schemes correct for a non-normal electron beam. In the program the three 
correction schemes are incorporated in separate subroutines so that they can be 
easily changed if better corrections are developed. 

1. Atomic Number Correction - The atomic number correction takes 
account of the fact that both the fraction of the total energy loss of an electron 
going into ionization of a particular shell (S) and the loss of ionization efficiency 
due to backscattering (R) are functions of atomic number. The atomic number 
correction is usually expressed as a ratio of backscattering to stopping power; 


[At. No.] u s' [R/S ] u 
[At. No.] g ft. [R/S3 s 


( 6 ) 


Of the 3 major corrections, the atomic number correction appears to be the least 
developed. Uncertainties appear to lie in the accuracy of the R values available, 
in the validity of the Bethe equation for stopping power and in the value of J to 
be used in the Bethe stopping pov/er. The correction according to Duncumb and 
Reed 4 provides a method for the correction that is easy to apply and appears to 
be achieving some acceptance. This formulation provides accuracies of the order 
of a few percent (Duncumb et al.). 14 


According to Duncumb-Reed, 4 


R = 


L R * c ‘ 


* 


to 
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for a complex sample where R l is the backs catter coefficient for each element, 
Cj is the mass fraction of the i th element and R i is a function of 1/U, the re- 
ciprocal of the overvoltage. S may be obtained by using an expression of the 
form 


where; 



$ A - Zj/A. In 


^I.166x 10 3 E ^ 


( 7 ) 


and E - (E 0 - E u )/2 for each element in the sample. E 0 is the operating voltage 
and E c is the excitation voltage of the x-ray line of interest. The mean ioniza- 
tion potential J v was determined empirically as a function of atomic number Z, 4 
The atomic number correction R/S can therefore be calculated for both unknown 
and standard. It should be noted that the correction is the same for all the ele- 
ments in the sample. 

2. Absorption Correction - The absorption correction for both the analyzed 
sample and standard is given by the absorption function f (X), where [Abs] u » 
f(X) u and [Abs] a = f (x) 8 . The absorption correction procedure of Philibert 1 
modified by Duncumb and Shields 2 is widely accepted as satisfactory, especially 
if f (X) has a value between 0.3 and 1.0 (Yakowitz and Heinrich). 8 

The absorption function for each element analyzed in a sample is given by: 

[Abs] u - f(X) u = I + h/(;i + X/c) [1 + h (1 + X/cr)] u ), 

( 8 ) 

[Abs] s = f(X) s = 1 +h/^l + X/cr) [1 + h (1 + X/o*)] s ) 


in which 


1. X = (p>/p) cscd where ji/p is the mass absorption coefficient of the 
target for the measured x-ray line; 

p/p = ' .Cj. p/pj < 9 ) 

i 
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and C i = wt.% of element i 

M/ Py ~ mass absorption coefficient of element i for the x-ray line of interest. 
6 = is the t alee -off angle of the emitted radiation. 


2, 


h r 1.2 


L 


a 



a i Zs 


( 10 ) 


where ai is the atomic concentration; A the atomic weight; and Z., the atomic 
number of element i in the multi-element sample. The atomic concentration, 
a. , is approximated by the weight fraction, C £ , within the program. 


3. 


, _ Cx 10 s 

\J •wt<i frjwM »i 


E, f n . E n 


( 11 ) 
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where E c is the excitation potential of the element of interest and E 0 is the 
electron beam operating voltage. Heinrich 3 proposed the values, C = 4.5 and 
n = 1,65. Recent studies by Duncumb, et al. 14 Show these values of C and n to 
be about optimum. The absorption function, f (X) , must be recalculated when 
another element is analyzed or when C A changes. This correction limits the 
general applicability of the calculation procedure because p/p values are not 
available for long wavelength L a radiation. 


3. Fluorescence Correction - The fluorescence correction for characteristic 
x-rays is required when one or more of the elements present in the specimen or 
standard has a characteristic emission line with energy greater than the critical 
excitation potential of the line being measured. In this case a fraction of the in- 
tensity observed, I f /X p , is due to fluorescence excitation rather than direct ex- 
citation. I f represents the intensity due to fluorescence and I p represents the 
primary intensity from the radiation of interest. The fluorescence correction 
for characteristic x-rays is given by: 


[f luor] u 


1 


[f luor] s 



( 12 ) 



6 


1 


where I f /I p is summed over all the elements which cause fluorescence. The 
secondary fluorescence correction procedure developed by Castaing 13 and 
modified by Reed s has received the widest acceptance and seems to be the most 
suitable for general usage. The modified version of Castaing's correction for- 
mula for K - K fluorescence involving L lines is given by: 


o 




= 0. 5P. .Ce yW(B) — 


r U B -l x1,67 


B\U a -1> 



ln(l+X) ln(l+Y) 
— X + y 


) 


(13) 


A is the element fluoresced 

B is the element whose characteristic radiation excites element 
A radiation 

C is the concentration of element B 

B 

r A is the absorption edge jump ratio of element A 

W(B) is the K or L shell fluorescence yield of element B 

A and B are the atomic weights of elements A and B . 

U B is the overvoltage ratio, E 0 /E c t for element B 

U A is the overvoltage ratio, E 0 /E c , for element A 

fjL* is the mass absorption coefficient of element A for radiation from 
B element b . 

/j,- is the mass absorption coefficient of the specimen for radiation 
from element B . 

The absorption parameters X and Y which appear in the above expression 
are given by: 


X= O a /> b ) csc&, Y = cr//z B 


(14) 


where / u, A is the mass absorption coefficient of the specimen for A radiation, 
o ■ is the electron mass absorption coefficient and # is the take-off angle of 
the emitted radiation. takes account of the extra factor required in K-L 
and L-K fluorescence, where j represents the excited line (K a , V L « > and i 
represents the excitation line (K a »K« »L a )* When only certain lines of R or L 
spectrum produce fluorescence excitation special precautions must be taken. 
The value of (r A — l)/r A for L radiation is equal to 0.75 (Reed). 5 
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The correction for fluorescence by the x-ray continuum has been developed 
by Henoc 15 and by Springer. 16 Generally it has been thought to be a small 
correction. Recently however, Brown, et al. 17 have shown that this correction 
is larger than often considered in the past. At this time our computer program 
does not include this correction. 


III. GENERAL DESCRIPTION OF THE PROGRAM 

A. Flow Diagram - The microprobe correction program is designed to run 
in either of 2 modes. Mode 1 allows one to compute the composition of the sam- 
ple of interest from measured data (I A and I (A). Mode 2 allows one to calcu- 
late the expected intensity ratios (I A /l (A)) from the known composition of the 
sample. The flow diagram for the correction program is given in Figure 7. 

1« Mode 1 •- The input data contains the parameters necessary for the 
calculation steps and the raw intensity data from both sample and standards. 

The raw intensity data is then compiled and the background, dead time and drift 
corrections are made. The mass absorption coefficients for all the elements 
analyzed and for all the elements which are present in both samples and stand- 
ards are then calculated using the method of Heinrich. 6 The corrections, ZAF, 

(Z is the atomic number correction, A is the absorption correction and F is the 
fluorescence correction) are then computed for the elements measured in each of 
the standards in turn (Equations f , 8 and 12). If any of the standards are pure 
elements (C(A) = 1.0) no secondary fluorescence correction, F, is made. 

The intensity data from the first sample data point is divided by the intensity 
from the corresponding element standard to obtain the measured intensity ratio, 
I A /I (A). This procedure is repeated for all the elements measured in the first 
sample data point. Then an initial estimate of the composition of all the elements 
measured in the sample, C° unknown , is made by assuming the linear relation 

Clknown = C A = C 1 ( A » measured C < A )' (15) 


If there are any elements in the sample which are not measured, their composi- 
tions are read into the program as input. If the sample contains oxygen, a cor- 
responding oxygen content is calculated from the assumed cation content of each 
element in the sample. The total amount of oxygen in the sample is then obtained 
to be used in the following calculations. The ZAF corrections are applied to the 
values of C° nknown and a calculated value of I A is obtained, (I A ) calc : 
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I 


* 


(U.u. = (ZAF) C u " nkn0W n 


(16) 


The calculated intensity ratio for each element of interest is then obtained, 
k CQlc = (I A /I (A)) calc . A new estimate of the compositionC A can then be deter- 
mined for each element from the equation: 


C' 

unknown 


(I ^/I (A)) measurec j 
(I a /1(A)) 

calculated 


X C 


o 

Unknown 


( 17 ) 


A test for the convergence of each element is made in the following fashion: 


if 


(V* (*)>«.! ciliated j 
(I ^/I (A)) measurec j 


<5x 10~ 3 


(18) 


sufficient convergence has occurred. If any of the elements do not converge, 
another iteration with new values of C^ nknown is made. Iterations are continued 
until convergence is achieved for all elements. 

When convergence is achieved, the calculated compositions are read out. 

The next point on the sample or another sample is selected and the calculations 
are repeated. It is only necessary to use new standards and re-calculate standard 
corrections when a new problem is introduced. 

2. Mode 2 — In this mode the compositions of both the standard and sample 
are input as data. One iteration of the correction procedure is made in order to 
calculate the intensity ratios, (I A /I (A )) calculated for each element in the sample. 

B. Structure of the Computer Program ♦ 

The program is written in the Fortran IV language according to specifications 

* for its use on IBM S/360 computers. It is built up in modular form using sub- 
routines each of which performs one task. A listing of the program and a glos- 
sary of terms used in the conversion program are given in the Appendices I and 

* II respectively. Detailed flow charts and descriptions of the most important sub- 
routines are found in Appendix HI. These can be used when the program is modi- 
fied or initially set up on another type of computer. The subroutines are described 
briefly in the following paragraphs. 
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MAIN - Organizes the input data, directs execution of the subroutines, calcu- 
lates cr (Equation 11), which is used in the absorption and fluorescence corrections, 
and coordinates the output of the final results. 

DATPTS - Associates for each of the measured elements, the standard back- 
ground and sample intensity values as well as the time for analysis. It also 
associates the intensity values with the spectrometers used and with the com- 
positions of the standards used for each element. DATPTS inputs the intensities 
measured by the probe to MAIN after calculating intensity in counts per second 
and is the one subroutine which is most dependent on the equipment which records 
the microprobe data (scalers, data translator, typewriter or tape readout). 
Therefore this is the routine which users will probably wish to tailor to their 
individual equipment. 

MODIFY — Computes the drift correction, if necessary, to be applied to the 
measured data. It also calculates dead time and background corrections for the 
measured standards and samples. 

ABSCO — Computes the mass absorption coefficients, according to Heinrich, 6 
for K a , Kp and L a radiation from all the elements present in both the standards 
and samples as absorbed by each of these elements. Some of these cannot be 
determined 6 and are calculated as zero. 

ATNCOR — Computes the atomic number correction for the standards and samples. 

INTER P - Interpolates linearly from the table for R given by Duncumb and Reed 4 
to obtain the value of the backscatter coefficient. The coefficient R is a function 
of both the atomic number of the element of interest, Z, and the reciprocal of the 
overvoltage, l/U. 

ABSCOR — Computes the absorption correction for any measured element in 
either sample or standard. 

FLUCOR — Computes the secondary fluorescence correction for K or L radia- 
tion produced by either K a , K^or L a radiation from any element present in either 
sample or standard. The total fluorescence correction (Equations 12-14) 

1+ E l < n * 


is obtained by the summation of the effect from the various elements in the 
sample. 


TEST - Determines whether secondary fluorescence by a sample matrix element 
can excite the element of interest. This is done by determining whether the K a , 
K /? and L a line produced by each element present in the sample is less than the 
K or L edge of the measured element. 


CHANGE - Calculates a new estimate for the composition of the unknown (Equa- 
tion 17) at the end of each iteration. 


ENDIT - Tests for convergence of each element in the sample at the end of each 
iteration (Equation 18). If convergence is achieved for all elements the composi- 
tion of the sample is printed out. 

OXYGEN - Calculates the amount of oxygen associated with each cation i if 
oxygen is present in the sample. The following equation is used to calculate the 
oxygen content: 


C for cation i = — — - — 
05f y A .OX. 


(19) 


k 


* 


Where OX. = N/M if the oxide has the stoichiometric formula A n O m . 
C . * wt. % of element i 
Ai == atomic weight of element i 

The total amount of oxygen present in the sample is given by 


C oxy (total) = 




where N is the number of measured cations in the sample and C o is the amount 
of oxygen associated with cations in the sample which are not measured. C o is 
given as input data. 


IV. INPUT TO THE PROGRAM 


The input data for the program includes all the necessary physical param- 
eters and compositional information for calculation of the corrections. If Mode 
1 is used the measured intensity data are also included as input. The method of 
input used for the program is a compromise between incorporating many of the 
parameters into the basic program, which uses up computer storage, and 
developing very flexible input which increases the time required for the analyst 
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to organize and run a problem. The input data are broken up into four files which 
are considered separately. 


A. Data File - This file contains the parameters which describe the indi- 
vidual problem and operating conditions of the electron probe. 

B. Element File - This file contains the parameters which are a function 
of a particular element. This file is already available and the elements 
which will be considered in the problem are selected by the analyst 
before the run is made. 

C. Composition File - This file contains the composition of the standards 
used for the analysis. The standards are usually available well in ad- 
vance of a particular run. The data contained in the composition file are 
often stored on cards or at the computer. 

D. Intensity Data File — This file contains the intensity readings from the 
probe. In many cases, punched tape or cards which are suitable for 
direct submission to the computer are obtained as output from the 
microprobe. 

The Element file and the Composition file are available beforehand, and in 
remote terminal usage individual file members can be selected and compiled 
into the data set for a problem just prior to executed time. The intensity 
readings from the probe are available or can be made available in suitable form. 
Therefore, the only input cards that need be punched, before a computer run can 
be made, are the DATA cards. Sets A, B, C are input in the MAIN program 
while subroutine DATPTS inputs set D. Appendix IV considers each of the 4 
types of input data separately and describes how the input must be set up to make 
a successful run. 

V. OUTPUT FROM THE PROGRAM 

We have selected the following output as most important. The element (E) 
file and selected parameters from the Data File are printed out first. Next the 
matrix of mass absorption coefficients are given, followed by the raw probe data 
as it was submitted to the computer. 

The correction ca„ illations for the standards are given next: (1) Atomic 
number correction for each standard, (2) Absorption correction for each element 
in each standard, (3) Fluorescence correction for each element in each standard. 
For the first point on each sample, the calculated element compositions are 
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printed for each iteration step and the ZAF corrections are printed for the final 
iteration step. For subsequent data points, only the composition values from the 
first and last iterations are written out. 

In Mode 2 problems, the calculated intensity ratios and the ZAF corrections 
are printed out for each point. In Appendix V the statements which control the 
output will be discussed as well as methods to change the output statements if 
desired. 
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APPENDIX I 
Program Listing 




MAIN PROGRAM 

S.N. 

^ ^ >H^< ^ ^ >§i ^ )i< J$e ^ ^ ^ ^ 1^ >50^e >S« 

100 

r. 



200 

c 



300 

c 

micrdprobe corrfct i on program march, i<?a9 

400 

c 



500 

c 



600 

c 


,!, 1*001.0$ TP IN 

700 

c 


DEPARTMENT OP METALLURGY fi MATERIALS SCIENCE 

BOO 

c 

\ 

LEHIGH UNIVERSITY 

900 

c 


R PTH L FH EM f p FNN 5 Y l. V A N I A 

1000 

c 



1100 

c 


P.A. COM 61 (.A 

1200 

c 


MATHEMATICS ANO COMPUTING BRANCH, COOP 64? 

1300 

c 


laboratory FOR THFORFTICAL studies 

1400 

c 



1500 

c 



1600 

c 


GODOARO SPACE FLIGHT CENTER 

1700 

c 


GREENBELT, MARYLAND 

1800 

c 



1900 

c 



2000 

C*^********^*****#**#^^^##****##******^**********#***********#****** 

2100 

3 

FORMAT ( • 

1ELEMENT* 3X ,A2 , 1 1 ( 8X , A? ) ) 

2200 

4 

FORMAT ( • 

ATOM. NO* ,1?F10*4) 

2300 

5 

FORMAT ( • 

ATOM.WT* ,12F10.4) 

2400 

6 

FORMAT ( ' 

K ALPHA 'tl2.Fl 0.4 

2500 

7. 

FORMAT 1 * 

K BETA * , l 2P10.4) 

2600 

8 

FORMAT! » 

K EDGE N12F10.4) 

2700 

9 

FORMAT ( • 

L ALPHA* ,1 2F10.4) 

2800 

10 

FORMAT ( * 

L3 EDGE* ,1?F10.4) 

2900 

ll 

FORMAT ( * 

CK • , 12F10.4) 

3000 

12 

FORMAT ( » 

NK *,12610.4) 

3100 

13 

FORMAT ( • 

CKL • ,12F10.4) 

3200 

14 

FORMAT ( • 

NKL * ,l?F10.4) 

3300 

15 

FORMAT ( 

* LI *,12610.4) 

3400 

16 

FORMAT ( ' 

CL1 *12610.4) 

3500 

17 

FORMAT ( » 

L2 • , 12F10.4) 

3600 

18 

FORMAT ( * 

CL? *,12610.4) 

3700 

19 

FORMAT ( * 

CLM *,12610.4) 

3800 

20 

FORMAT ( * 

Ml *,12610.4) 

3900 

21 

FORMAT ( • 

CMl • ,12610.4) 

4000 

22 

FORMAT ( ' 

M2 • 12F10.4) 

4100 

23 

FORMAT ( ' 

CM2 * ,12610.4) 

4200 

24 

FORMAT ( ' 

M3 *,12610.4) 

4300 

25 

FORMAT ( • 

CM3 ' 12F10 .4 ) 

4400 

26 

FORMAT ( • 

M4 *126)0.4) 

4500 

27 

FORMAT ( • 

M5 ' 12F10.4) * 

4600 

28 

FORMAT ( * 

CMN *12610.4) 

4700 

29 

FORMAT ( • 

N1 '12610.4) 

4800 

30 

FORMAT ( ' 

VC(Z,K) *12610.4) 

4900 

31 

FORMAT ( ' 

VC ( Z ,L ) ' 12610.4 ) 

5000 

32 

FORMAT ( • 

J1IZ) *12(18, 2X)) 

5100 

33 

FORMAT ( • 

RK(Z) *12610.4) 

5200 

34 

FORMAT ( ' 

OXY(Z) *12610.4) 

5300 

35 

FORMAT ( • 

OMEGA-K* 12610.4) 

5400 

36 

FORMAT ( ' 

OMEGA-L* 12610.4) 

5500 

37 

FORMAT ( • 

LINE • 3X,A2 ,11 ( RX,A? ) ) 

5600 

38 

FORMAT ( ' 

VO '126)0.4) 

5700 

39 

FORMAT ( • 

SPEC. NO* , 12(17, 3X) ) 

5800 

40 

FORMAT ( 20A4 ) 

5900 

41 

FORMAT ( 1 4 , 10( ?;X,A2 ) ) 

6000 
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p 


4? 

43 

44 

45 


46 

47 
4fl 
49 


50 

51 

67 


705 

706 

707 

710 

711 
713 

?14 

715 

716 

717 

718 

719 


720 


C 


FORMAT < 4X • 1 PF4.2 ) 

FORMAT ( * THFRF ARP ' ,13,' FI.PMFMT 5 |M TH I 5 SAMPI.F !* 5X • A2 » 9 ( 6X , A2 


NUMBER CORRECT ION» ' » P13.5 ) 

FOR * , A2 » ' CHI*' ,F13.5,7X, 


1 )) 

FORMAT! • OATOMJ C 
FORMAT! 

1 / ‘OAR SORPTION CORRECT I ON 

2 'FICHU-' | F l ? • 5 ) 

FORMAT! » WEIGHT PERCENTAGES /100 FOP TMF FI FMFMTS * ' P7 .4 , 9FR.4) 

FORMAT!////, 'O' ,20AM 

FORMAT !////• OF I NAI. CORRECTIONS'//) 


FORMAT! 


1 

2 
2 
2 
2 
2 


THFTA 

TAII 

N 

NM 

NR OS 
MOOF 
NPROR 


• , FI 0.4 / 
' 3E13.5/ 
17/ 

17/ 

17/ 

17/ 

m 


FORMAT! • 1' ) 

FORMAT! 'OELEMENT ',A?,» IS MOT ACCOUNTED FOR— NOT INPUT, MOT 
1 MEASURED, MOT OXYGEN') 

FORMAT (IX, A 2, IX, 

1 I 2 ,2X ,F5. 1 ,2 ( IX ,F6. 3 ) , IX, 1 4, 1 X »F5 .2 » IX » F6.4 , l X * F5.3 » 3 (1 X 

1,F6.3) ,1X,F7.374X,F4. 2,1X,F6.3, IX, F6. 2, IX, F4. 2, 1 X, F6.4, 1 X, F6. 3, 5! 1 


2X,F6.2)/4X,F5U ,2(IX,F6.3,1X,F4.1) ,2! IX,F6.3> , IX , F5.2 , IX , F5.3 , IX , 
3F6.5) 

FORMAT! *0 ATOMIC NUMBER CORRECTION | ANC ) — ANC a ATN/ATO= ( RRU/SR 

1U)*(SRS/RRS) -( • , A2 , • ) • I 

FORMAT ( ' ' » 6 X , » A R S =' ,F10.5/7X, ' ARM s' , FIO .5/7X , ' ARO «',F10.5) 

FORMAT! '0 ARSORPT I ON CORRECTION! ARS ) — ARS-ARN/ARO ! • ,A2, 

l»l ' ) 

FORMAT (7X , ' ATMs • Fl 0. 5/7X, ' ATO=* FJ O. 5/7X , 'RRU® ' F I 4. 5/7X ' SRU= ' F14. 5 
1 /7X , ' RBS = ' F14.5/7X ' SRS a 'El4. 5/// / ) 

FORMAT !7X,' AMC®' ,F10.5) 

FORMAT! '0 FLUORESCENCE CORRECTION! FLO) — FLU=F!,N/FI,0 C,A 

12, • )- • ) 

FORMAT!' ' ,6 X, 'FLU = » ,F10.5/7X, • FLM s' , FIO .5/7X , ' FLO =' , E10.5 ) 
FORMAT!/////* ELEMENT CHART — WEIGHT PFRCFMT ' / )X,64( •-• ) ) 
FORMAT!' DATA PT NO ITER ' , 3X, A2 „ 10 ( 7X , A3 )) 


FORMAT ( 1 5 , 5X , I 5 , 3X , 1 1 ( 2PF10.3) ) 

FORMAT! 15, 6X,A4,3X,1U2PF10. 35 1 

FORMAT! 'OMODE 2 PROBLEM: INPUT-COMPOSITION R AT I OS. . . OUTPUT-! NTENS I 
1TY RATIOS'/' PROBLEM#' ,15/' ELEMENT.. . . . . INTENSITY RATIO-SAMPLE TO 
3STAN0ARD-ATN*ARN*FLN*C(U)/ATn*AP0>!'Fin*C(S)'/(4X,A4,10X,Fl?.5) ) 
FORMAT! 'OIN STANDARD COMBINATION #• I5,2X,A2. • USFS STANDARD #*13) 
DIMENSION nx Y ( inn) ,S10I9,10> ,AK(9), 

1 AK 1(10) ,UC (JO) ,SC ( 10) , SCON (9, 10) ,ATDS(10) ,AR0S(9,10 ) , FLOS (9, 10) , 
l RRS(IO) tSRS(lO) 

INTFGER*4 STD,ll,RSTO 

RFAL*4 NK,WKL,L1,L?,L3,M) ,M2 ,M3,M4 t M5,N) ,K ALPHA, KRETA,KFDGF, I. ALPHA 


NAMFLIST /DATA/ THFTA , TAU, STO ,N , NFII,MOnF , VO, L I ME, I SPEC , NPROR , 

l nstcom,mm,nrgs,rstd,m 

COMMON/EV/VC (100,2 ) /JJJ/Jl! 100 ) /RVAI./RK ( 100 ) /WVAL/OMEGA ! 100,2) 
A /F9/ ARCO! 50, 50, 3), SIG, THETA 


X 

Y 

X 

X 

1 

2 

3 

4 

5 
X 
X 


COMMON / RACK/ RGINT! 9 ,9 ) ,RKGR ( 9 ) , BSTO! 9 ) ,NRGS , RTITI.F! 20,9 ) 
/STAND/ ESTO! ?0,10) ,ST0INT(9,9,).0) ,NESTO! 10), 

CS( 10,10) ,STINTD(9,10,P) ,C0MPST(9,10) 

/SCRIPT/ N ,NU( 8 ,3 ) ,NM »M , INDTC 

/INFO/ E ( 50 ) ,Z ! 50 ) , A ( 50 ) ,KALPHA ( 50 ) ,KRFTA ( 50 ) ,KFOGF( 50 ) 

L ALPHA ( 50 ) ,L3 ( 50 ) ,CK( 50 ) ,NK ( 50 ) , CKL ( 50 ) ,NKL ( 50 ) , 
L1!50),CU (50), 1.2(50), CL 2! 5n ) ,CI M( 50 ) ,M) ( 50 ) , 
CM1!30),M2(50),CM2(50),M3!50),CM3(50),M4(50), 

M5! 50 ) ,CMN ( 50 ) , N1 ( 50 ) 

/ARCS/ LINF(IO) ,VO,MP.I , L |. • L LL , , NCOHMT 
/CORREC/OR I FT! 9, 1 0 ,B ) , TAU(3), JSPFC(9) 

/IINKWN/ HINT! 75,9 tfl) ,NEH,C( 10) 
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6100 
6200 
6300 
6400 
6500 
6600 
6700 
6ROO 
6900 
7000 
7100 
7200 
7300 
7400 
7500 
7600 
7700 
7B00 
7900 
BOOO 
R 100 
8200 
8300 
R400 
8500 
8600 
8700 
8800 
8900 
9000 
9100 
9200 
9300 
9400 
9500 
9600 
9700 
9800 
9900 
10000 
10100 
10200 
10300 
10400 
10500 
10600 
10700 
10800 
10900 
11 000 
11 100 
11200 
11300 
11400 
11500 
11600 
11700 
11800 
11900 
12000 
12100 
12200 
12300 
12400 
12500 
12600 
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n 


X /XXX/ CHI 12700 

a /inoiit/ iout.jn 121100 

RFAL*4 ST I TLF ( ?n, 10) ,IIT!TI.F(?0) 12900 

SIGMA (Q,W)» 4 , 5P5/ ( 0**1 ,65-W#*J .AS) 13000 

imiT»A moo 

IN**) 13200 

DATA LAST .L »n, SUM/ * LAST * . < L ' • * n* * ' SUM » / 13300 

100 RFAD ( IN f DATA »FNn*l 0000 ) 13400 

NIMlti)«l 13*500 

INDlC*l 13000 

NPl *N+ J 13700 

IS*0 13B00 

00 i?5 1*1, M 13900 

READ! IN » 67 1 RND* 10000 ) FID U ,VC! I , l ) .VC ( 1 ,2 I ..II (!) ,RK l ] ) , 14000 

1 nXY I I ) fDMFGA (1,1) 14100 

1 ,KA1PHA ( I ) ,KBFTA ( I ) ,1A1PHA( I ) ,CK( I ) ,NK( I ) ,KFDGF ( I ) ,f.KU I ) ,N«U I ) • 14200 

miihcuii ) ti.?( n tCL?u)»L3U)»ci.M( i > ,mi m,cMi 14300 

3 M3intCM3(I),M4U)*M<5m,CMN{I),Nim,nMFGA(I,?) 14400 

7. ( I ) *FLOAT I NO) 14500 

125 CONTINUE 14600 

MP*M 14700 

1*1 14B00 

0*1? 14900 

130 IF(M-U) .135,135,160 15000 

135 U*M 15100 

160 WRITE ( I OUT ,3) <F(I),1*1,U) 15200 

WRITE! I OUT ,4} < Z 1 1 ) , J *1 ,0) 15300 

WRITE 1 1 OUT ,5) (Al I ) ,1*1,11) 15400 

WRITE! 1 OUT, A ) ( K ALPHA ( I ) , I *1 »U ) 15500 

WRITE ! IOUT ,7 ) I KBETA ( I ) , I = 1 ,tl ) 15600 

WRITE! IOUT,B) (HEDGE! 11,1*1,11) 15700 

WRITE ! IOUT ,9) UA1PHA ( l ) ,1*1. ,U) 15R00 

WRITF! i0WT,10ML3U),I*L,U) 15900 

WRITE! IOUT, 11) (CK( I ) ,I*1,U> 16000 

WRITE! I OUT , 1? ) ( NK ( J),I#L,U) 16100 

WRITE! IOUT ,13 ) ( CK1 ! I ) ,1*1,11) 16200 

WRITE! IOUT ,14) (NKK I ) ,1*1,0) 16300 

WRITE! I OUT ,15)! 1.1 ( 7 ) , l*L,U) 16400 

WRITE! I0UT,16)(CU!I ) »I=L,U) 16500 

WRITE! I OUT, 17) <l?.( I ) ,1*1,11) 16600 

WRITE! |OUT,lfl)(Ct?( I ) ,1*1,11) 16700 

WRITE! I OUT, 19) (CUM! I) ,1=1,11) 16R00 

WRITE! IOUT ,20 ) ( Ml ( I ) , I *1. ,11 ) 16900 

WRITE! IOUT, 21 ) !f.Ml! I ) ,1=1, U) 17000 

WRITE! IOUT,??) (M?< 1),I*1,U) 17100 

WRITE! IOUT, ?3) ! CM2! I ) ,I=1,U) 17200 

WRITE! I OUT, 24HM3! !),!=(., ID 17300 

WRITE! IOUT ,75) I CM3( I ) ,1*1,0) 17400 

WRITE! I OUT, 26) ( M4( 1 ), I*l,!l) 17500 

WRITE! IOUT, 27) (M5! I ) ,1*1,11) 17600 

WRITF ( IOUT , 2R ) ( CMM ( I ) ,1*1, U) 17700 

WRITE! IOUT ,29 ) ! NH I ) , I *1 ,U) 17R00 

WRITE! IOUT, 30) ( VC ! I , 1 ) ,1*1,11) 17900 

WRITE! IOUT, 31) ! VC ( I ,? ) , I =1. ,11) 1B000 

WRITE! I OUT, 3?) ! Jl ( I) , 1*1, ID 1R100 

WRITE! IOUT, 33) ( R K ! I ) ,1*1,11) IR200 

WRITE! IOUT, 34) (OXY(I), 1*1,11) IR300 

WRITE! IOUY, 35) (OMEGA! 1 ,1) , 1*1,11) 1R400 

WRITE! IOUT, 36) (OMFGA (1 ,2 ) , 1=1,11) 1R500 

WRITE! inUT,37MlINF( I ) ,I=1,U) IR600 

WRITE! IOUT , 39 ) ( I SPEC! I ) , 1 = 1 ,N ) 1R700 

IF(U.GE.M) GO TO 175 l«800 

1*U+1 1 B900 

U*U+1? 19000 

WRITE! 6,41) U 19100 

WRITE (6 ,41) M 19200 

GO TO 130 19300 
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175 CONTINUE 19400 

WRITE! I OUT ,3R) VO 19S0Q 

WRITE! 1 01JT *49) TMFTA , ! TAUII) , 1 *1 ,3) .M f MM t MAr,S ,MflDF 19600 

l ,NPROB 19700 

THETA»THETA*3.lAl«j9?6/lA0.0 19R00 

10X1=0 19900 

DO 200 1*1 ,NFU 20000 

200 IF (E(I).EQ.n) on TO 300 ?0100 

GO TO 400 20200 

300 10X1*1 20300 

400 CONTINUE 20400 

REAO( IN, 40 ) (UTI'TLF! I ) , 1*1 ,70 ) 20500 

READ! 5,2 ) ! G ( I ) , I * 1 , NfilJ ) 20600 

2 FORMAT !8F l 0.7) 20700 

DO 500 I *1 ,NFU 20R00 

500 C( I )* .01 AC (IT 20900 

IF! I0X1.GT.0) COXY* C< 10X1) 21000 

DO 600 K«1,NM 21100 

READ! IN, 40) f STITLF! I ,K) ,1*1,20) 21200 

READ! IN, 41) NEST, tESTD! 1 ,K ) , I»l ,NF$T) 21300 

READ! IN ,2 ) ICS! I ) , 1*1 ,NFST ) 21400 

DO 550 I*S,N6ST 21500 

550 CS!I,K)*.01*CS!I,K) 21600 

NESTD(K )»NEST 21700 

600 CONTINUE 21R00 

IF! NBGS.EQ.O ) GO TO 900 21900 

DO BOO K*1 ,NRGS 22000 

READ! IN, 40) (RTITI.fi!!, K), 1 = 1 , 20) 22100 

ROO CONTINUE 22200 

900 CONTINUE 22300 

IF(NPl.GT.MEU) GO TO 920 22400 

DO 910 I*NPl,NFII 22500 

IF(E( I). EO.O) GO TO 910 27600 

IF(C! I ).6E.l.E-5) GO TO 910 22700 

WRITE! I OUT ,51 ) F(I) 22R00 

STOP 22900 

910 CONTINUE 23000 

920 DO 930 1*1, M 23100 

1 JK=1 23700 

I F( L JNF! I ) . EO.L ) UK = 2 7^300 

930 LINE! I ) * IJK 3^400 

CALL ABSCO 23ROO 

IF(MODE.EO.l) CALL DA TRY 5 73600 

NC0UNT*1 23700 

DO 4000 K»t ,NM 23800 

NEST=NFSTO(K) 23900 

CALL ATNC0R(RBS(K),SRS(K),NEST,CS(1,K),ATDS(K),EST0(1,K) ) 24000 

WRITE! I OUT ,47 ) ( ST ITLE ( I .K ) , I =1 ♦ 20 ) 24100 

WRITE! 10IIT, 43) NEST , ( ESTD( 1 ,K ) , !=1 »NFST ) 74200 

WRITE! I OUT, 46 ) (CS( I ,K) ,1 = 1 ,NfiST) 24300 

WRITE ( I OUT, 44) ATDS(K) 24400 

DO 3500 1=1, M 24500 

DO 3500 J*1,NFST 74600 

IF (EUI.Nfi. ESTDU,K) ) GO TO 3500 24700 

LLL=Z( I ) 24R00 

DO 3200 J J= 1 ,NSTCOM 24900 

3200 IF! StO! I ,JJ ) .EO.K ) GO TO 3300 25000 

GO TO 3500 25100 

3300 LL*L INE { 1 ) 25200 

IF! LLL.GE. 37 ) LL=? 25300 

SIG=SIGMA(VO,VC( I , L I. ) ) 25400 

SCON! I,K)=CS(J,K) 25500 

CALL ABSCOR (NEST, CS ( l,K ), ARDS! I,K),E.STD(l,K) ) 25600 

WRI TE! IOUT, 45) F. ( I ) ,CHI ,ARDS ( I , K ) 25700 

. FLOS! I ,K )=1 • 25800 

IF! NEST.GT. 1 ) 25900 

1 CALL FLUC0R!CS(1,K),FLDS(1,K),ESTD(1,K),NEST) 26000 
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3500 CONTINUE 
4000 CONTINUE 

DO 7000 KKKK=1,NSTC0M 
WRITE! IOUT *50 ) 

DO 7000 I GROHP= I y INDJC 
NUP-NU! I GROUP » 1 J 
DO 7000 NC= 1 »NUP 
GO TO (4050,4100) , MODE 
4050 NCOUNT-O 

GO TO 4200 
4100 NCOUNT-1 
4200 DO 5050 1 = 1, N 
K«STD( I,KKKK) 

SC! I )-SCQN( I ,K) 

IF(NC.EO.l) WRITE! IOUT, 720) KKKK,E(I),K 

I F ( MODE .EQ .2 ) GO TO 5050 

DR-DRIFT! I ,K , I GROUP ) *NC 

ST INT-ST INTO! I ,K , I GROUP ) +DR 

AM I I-UINT ( NC, I , [GROUP ) /STINT 

C( I ) *AK ( I )*COMPST ( I ,K ) 

5050 CONTINUE 

5075 IF (IOXl.NE.O) CALL OXYGEN! OXY,C , A , N,lOXl,COXY) 
SUMM-0.0 
DO 5150 J=1,NEU 
UC( J)-C(J) 

5150 SUMM-SUMM+UC! J) ' 

INDEX-0 

IF(NC.GT.i) GO TO 5?00 
WRITE! IOUT, 715) 

WRITE! IOUT ,716 ) HI! I ) ,1=1 ,NEU) ,SUM 
5200 CONTINUE 

.WRITE! IOUT ,717) NC, INDEX, (UC( I ),I=1,NEU> ,SUMM 
5250 DO 5500 INDEX-1,10 
472 CONTINUE 

CALL ATNCOR < RRU,SRU,MEU, UC , ATN , F ) 

DO 5400 1 = 1 , N 

K=STD( I ,KKKK ) 

5325 LLL-ZU) 

LL=L INE! I ) 

IF!LLL.GE.37) Ll.=2 

sigbsigma! vo,vr,( i ,i.u ) 

CALL ARSCOR ( MEU ,IIC , ARM ,E ) 

CALL F LUC OR ( I !C, , FLN , E , MFIi ) 

ANC=ATN/ATDS(K) 

ABS-ARN/ARDS! I »K ) 

FLU-FLN/FLDS! I ,K) 

AAA=ANC#ARS*FLU 

AK1! I ) = AAA*UC( I )/SC( I ) 

IF(NCOUNT.EO.O) GO TO 5400 
WRITE! IOUT, 713) E(I) 

WRITE! I OUT, 714) FLU, FLN, FLOS! I ,K ) 

WRITE! IOUT, 707) F.(I) 

WRITE! IOUT, 706) ABS , ARN , ARDS ( I , K ) 

WRITE! IOUT, 705) Ed ) 

WRITE! IOUT, 711 ) ANC 

WRITE! IOUT, 710) ATM , ATDS ( K ) ,RRU, SRU, RRS ( K ) » SRS ( K ) 
5400 CONTINUE 

I F ( MODE • EQ.2 ) GO TO 5550 
IF! NC0UNT.ro. L) GO TO 5575 
CALL CHANGE ( N, AK1 , AK $UG ) 

IF (IOXl.NE.O) CALL OXYGEN ( OXY , f 1C , A , M, IOXJ ,COXY) 

IF(NC.NE.l) GO TO 5475 

SUMM«0. 

00 5450 1 = 1, NEIJ 
5450 SUMM=SUMM+UC ( I ) 

WRITE (6,717) NC , INDEX » ( UC ( I ) , 1=1 ,NEU ) , SUMM 


26100 

26200 

26300 

26400 

26500 

26600 

26700 

26800 

26900 

27000 

27100 

27200 

27300 

27400 

27500 

27600 

27700 

27R00 

27900 

28000 

28100 

28200 

28300 

28400 

28500 

28600 

28700 

28800 

28900 

29000 

29100 

29200 

29300 

29400 

29500 

29600 

29700 

29ROO 

29900 

30000 

30100 

30200 

30300 

30400 

30500 

30600 

30700 

30800 

30900 

31000 

31100 

31200 

31300 

31400 

31500 

31600 

31700 

31800 

31900 

32000 

32100 

32200 

32300 

32400 

32500 

32600 


t 
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5475 CAL.L ENQ1T(N, A*<1 ,AK,f.55?P) 32700 

5500 CONTINUE 32800 

5525 IF(NC.NE.l) GO TO 5600 32900 

NC0UNT=1 33000 

WRITE! IOUT, 4R) 33100 

GO TO 472 " 33200 

5550 WRITE! 6,719) NPROR , ( F ( I ) , AK 1 ( I ) , I = ) ,N ) 33300 

GO TO 7000 . 33400 

5575 IF! NC.NE. 1 ) GO TO 5600 33500 

WRITE! IOUT, 715) 33600 

WRITE! IOUT, 716) ! E! I) , I *1 ,NEU) ♦ SUM 33700 

5600 IF (IOX1.NE.O) CALL OXYGEN! OXY,UC , A , N, 10X1, COXY) 33800 

SUMM«0.0 33900 

DO 5625 1=1, NEU 34000 

5625 SUMM=SUMM+UC ( I ) 34100 

INDEX=LAST 34200 

5650 WRITE! IOUT, 718) NC , .’NDEX, < UC ( I ) , I = 1 ,NF.I) ) ,SUMM 34300 

7000 CONTINUE 34400 

GO TO 100 34500 

10000 RETURN 34600 

END 

SUBROUTINE ABSCO 100 

REAL*4 KALPHA »KRETA ,KEDGF ,1. ALPHA , LAMBDA 200 

REAL**' NK,NKL,L1 ,L?,L3,M).,M?,M?,M4,M5,N1 300 

COMMON /SCRIPT/ N,NIJ! 8 , 3 ) ,NM,M 400 

X /INFO/ E ( 50 ) ,Z I 50 ) * A ( 50 ) , K A LPHA ! 50 ) ,KRET A ( 60 ) , K FDGF ! 50 ) ♦ 500 

1 LALPHA ( 50 ) ,L3( 50.) ,CK ( 50 ) ,NK ( 50 ) ,CK|. ( 50 ) «NKL ! 50 ) , 600 

2 L l ( 50 ) ,CI. 1 ( 50 ) » L? ( 50 ) *CL? ! 50 ) »CLM { 60 ) ,M) ( 50 ) , 700 

3 CM1 ( 50 ) ,M? ! 50 ) , CM? ! 50 ) , M? ( 50 ) , f,M? ! 50 ) , M4( 50 ) , 800 

4 M5< 50) ,CMM( 50) ,M) 1 50 ) 900 

X /F9/ ARCO! 50,60,?) /IMDIIT/ I HI IT 1000 

DIMENSION LAMBDA ( 50 ) 1100 

C THIS SUBROUTINE COMPUTES AND ASSOCIATES THE MASS ARSORPTION COEFFI 1200 

ECU XI ,Y1 ,Z1 )=Xl*Yl**Zl 1300 

WRITE! IOUT. 70) ) l Ano 

DO 80). JUMP® 1 ,3 1500 

GO TO (841, 844,84?), JUMP 1600 

841 DO 851 J = 1,M 1700 

851 LAMBDA! J)=KALPHA(J) l p0 ° 

GO TO 814 l 9no 

84? DO 85? -J=1,M 2000 

85? LAMBDA! J)=KRFTA(J) 2100 

GO TO 814 2200 

844 DO 854 J=1 ,M 2300 

854 LAMBDA ! J ) =L ALPHA ( J ) 2400 

814 DD R01 1=1, M 2500 

DO 801 J=1,M 2600 

ABCO! I , J » JUMP ) »0 .0 2700 

IF! LAMBDA ( J ). EO.O. ) GO TP 801 2800 

IF ( LAMBDA (J).LT.KFDGE(I ) ) GO TO 1 2900 

IF! LI ( 1 ) .EO.O.O ) LI ( I ) = 50,0 3000 

1 F (LAMBDA (J).LT.Ll(I) ) GOTO? 3100 

IF (LAMBDA! J).LT.L?( 1 ) ) GO TO 3 3200 

IFILAMBDA! J) .LT.L3! I ) ) GO TO 4 3300 

IF! LAMBDA (J).LT.Ml! I) ) GO TO 5 3400 

IF (LAMBDA! J) .LT.M2! I ) ) GO TO 6 3500 

IF! LAMBDA ( J ) .LT.M3! I ) ) GO TO 7 3600 

IF!LAMBDA(J)*LT»M4(I) ) GO TO 8 3700 

IF! LAMBDA ( J ) • LT.M5! I ) ) GO TO 9 3800 

IF(LAMBDA(J) •LF.Nl! I ) ) GO TO 10 3900 

GO TO (1000, 1100, 1000), JUMP 4000 

1000 WRITE! IOUT, I?) LAMBDA ( J ) ,F ( J ) ,N1 ( I ). F:( 1 ) 4100 

GO TO 5 4200 
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lino write( rniJT.m lamrda! j) ,E(jj ,nk n ,pi i ) 4300 

go rn 5 44no 

1 ABCO! I ,J,JUMP)sCK( I IHAMRDA! J)**NK( J ) 4500 

GO TO POl 4600 

? ABCfM I ,J ,JUMP)=CKI.( I)*LAMRnA( J)**NKI.U) 4700 

GO TO 801 4P00 

3 ABGnU,J f JUMP)sr.Ll( I )*LAMROA( .J)«*NKI.U) 4900 

GO TO POl 5000 

4 ABCOH ,J , JUMP > =CL? ( I )#LAMRDA( J )»KMK!,( I) 5100 

GO TO POl 5200 

5 ABCD! I ,J,JUMP)=CI.M( I)*LAMRDA(J )•!<*?. 60 5300 

60 TO POl 5400 

6 ABCn(I,J,JUMP)=CMl( I )*LAMRDA( J)**?.60 5500 

GO TO POl 5600 

7 ABCOU r J,JIJMP)=r,M?( nnAMROAUI) W.60 5700 

GO TO POl 5P00 

fl ARCO! I «J , JUMP ) =CM3 ( I ) *LAMRD A ( ,1 ) **? . 60 5900 

GO TO POl 6000 

9 WR I TE ( 6 » 13 ) LAMRnA ( J ) 6100 

GO TO POl 6?00 

10 ABCO( I,J,JUMP)=CMNU)*LAMROA( J)**?.?? 6300 

POl CONTINUE 6400 

WRITE! I OUT *701) 6500 

00 901 JUMPS] f 3 6600 

GO TO ( 1200,1300,1400) , JUMP 6700 

1200 WRITE! I OUT, 70?) 6P00 

GO TO 1500 6900 

1300 WRITE! I OUT \ 706) 7000 

GO TO 1500 7100 

1400 WRITE! IOUT, 705) 7?00 

1500 L= 1 7300 

MI=MIMO( I 4,M) 7400 

1550 WRITF! IOUT, 703) ( F ( I ) , I s|. ,M I ) 7500 

00 ) 600 1=) ,M 7600 

1600 WRITE! IOUT, 704) F ( I ) , ( ARGOd , J , JUMP ) . J = |. ,M t ) 7700 

IF (MI.GF.M) GO Tn 901 7R00 

L =M I + 1 7900 

MI =M I NO ( M I + 1 4 , M ) POOO 

GO TO 1550 R ] 00 

901 CnMTIMUF R? 00 

701 FORMAT! ' lCnRRFCTIOM VAl.UF.S (FIRST nATA POINT , I.AST ITFRATTOM R300 

1 ) ‘/IX, 57! ) ) R400 

702 FORMAT!' MASS ARSORPTIOM COFFFTC IFMTS (FMITTFR AT TOP) — K R500 

1 ALPHA RADIATION' ) R600 

703 FORMAT! 14(7X, A?) ) R700 

704 FORMAT! 1X,A?,?X,14(F7.I ,?X) ) RBOO 

705 FORMAT! ' MASS ARSORPTIOM GOFFF T C I FMTS ( EM I TTFR AT TOP) — K R900 

1BFTA RADIATION' ) 9000 

706 FORMAT!' MASS ARSORPTIOM COFFFICIFNTS (FMITTFR AT TOP) — I. 9100 

1 ALPHA RADIATION') 9200 

12 FORMAT!' LAMRDA-K ( • , F6 . 3 , ' ) OF 'A4,» IS GREATER THAN N](',F6.3, 9300 

1 • ) OF ' A4 ) 9400 

13 FORMAT!' LAMRDA-L ( • *F6.3, ' ) OF *A4.» IS GREATER THAN NU'»F6.3t 9500 

1 ' ) OF *44) 9600 

END 9700 


SUBROUTINE ABSCOR (NEL,C,AB,EL) 100 

RE AL’f‘4 6L< 10) ,C( 10) 200 

COMMON/XXX/X 300 

1 /INFO/ E ( 50 ) »Z ( 50 ) » A ( 50 ) 400 

3 /F9/ ARCO! 50,50,3) ,SIG, THETA 500 

4 /ARGS/ LINE! 10) ,VO,M,I ,KK 600 

AAVE=0.0 700 

ZAVEsO.O BOO 

X=0.0 900 
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DO 2 .1=1, NEL 
on ? k=i,m 

IF(F(K).NE.EL(.l) > GO TO 2 
AAVF*AAVE+C( J)*A(K) 

ZAVB*7.AVE+CI.I)*7(K) 

X = X+C<.l)*ARCO(K,I,KK) 

2 CONTINUE 

H*1.2*AAVE/ZAVE#«2 

X=X/SIN(THFTA) 

ar=( i.o+H)/m.n+x/sifi)#i i.o+h*i i .o+x/sig) ) > 

RETURN 

.END 


iono 

1100 

1200 

1300 

1400 

1500 

1600 

1700 

1R00 

1900 

2000 

2100 


SUBROUTINE ATNCOR (RB,SB,NEL > C,AT,£L) 100 

DIMFNS ION EL ( 10 ) ,C ( 1 0 ) 200 

COMMON/EV/VCl 100,2)/. KU/JI 100) 300 

1 /INFO/ F( 50), l { 50), A (50) 400 

2 / A RGS / L I N F ( 1 0 ) , VO ♦ M 500 

RR=0.0 600 

SR =0,0 700 

DO 1000 1=1, M ROO 

DO 500 K = 1 ,NFI_ 900 

500 JFUM I J.EO.EUK) ) GO TO 600 1000 

GO TO 1000 1100 

600 L*Z< I > 1200 

LL*LINEI I ) 1300 

I F ( L ,GF,37 ) LL=? 1400 

V=.5#IV0+VC( I »LL) ) 1500 

RECU=VC( I ,LL)/VO 1600 

CALL INTERP(RECU.L»RAR) 1700 

RB=RR+C(K)*RAR 

S = Z( I )/A( I )*AI. "IG( ( ) ,166F3*'//,)| III) 1900 

SR=SR+C(K)*S 2000 

1000 CONTINUE 2100 

AT=RR/SR 2200 

RETURN 2300 

END 2400 

SUBROUTINE CHANGE (N.AK1 ,AK,UC) 1 0 0 

DIMENSION AKl ( 10) ,AK ( I(ir,UC'UOT 200 

C THIS SUBROUTINE GIVES NEW VALUFS FOR THE WF I GHT PERCENTS 300 

DO 200 1=1 ,N 400 

J=I 500 

100 RAT=AK ( I ) /AK1 ( I ) 600 

UC< J)=RAT*UC(.I) 700 

200 CONTINUE ROO 

RETURN 900 

END 1000 


SUBROUTINE DATPTS 100 

TNTFGER*4 BSTD i SAME 200 

RE AL*4 TBG ( 9 ,3 ) ,TSTD( 9,10,3) , X INT ( 3 ) 300 

COMMON/SS/S ( 3 ) 400 

COMMON /BACK/ BGINT(9,9) ,BKGR ( 9 ) , RSTD( 9 ) ,NRGS 500 

X /SCRIPT/ N,NU(R,?) ,NM,M, IMOIC 600 

X /STAND/ ESTDI 10,10) , STDtNT (9,9,10) ,NFSTD ( 1 0 ) 700 

X /UNKWN/ UIMTI 75,9,P ) 800 

X /INOUT/ I OUT , IN 900 

X /INFO/ E (50) 1000 

WRITE! IOUT ,701 ) 1100 

DO 100 K= 1 ,NM 1200 

DO 50 L = 1 ,9 1300 

DO 40 J*1 ,3 1400 

TSTDIL,K,J>=0, 1500 
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40 CONTINUE 

DO 50 ,1=1, N 
50 STOINT ( L ,vl ,K ) =0. 

100 CONTINUE 

JE(NBGS.EO.O) GO TP ?75 
DO 250 K=1,NBGS 
DO 200 J=l,3 
TBG( K , J ) =0 . 

200 CONTINUE 

DO 250 J = i *N 
BGI NT ( J,K)*0. 

250 CONTINUE 
275 CONTINUE 

DO 2R0 L=1 *ft 
DO 2B0 1=1,3 
2B0 NU(L , I ) =0 

DO 20000 I CT= 1 ,2000 

3000 READ( IN , 1 ) TI MF ,X 1NT ,.SAMF,NIJMRDG, IMDIC ,NOST, 1 1 INK , I AST 

GO TO ( 3025,3050,3075) ,NHMROG 
3025 I T=1 

I NDMAX=M I NO ( N , 3 ) 

GO TO 14000,5000,(5000) , SAME 
3050 JT=2 

NUMRDG=4 

INDMAX=MIN0(N,6) 

GO TO (4000,5000,6000) , SAME 
3075 JT-3 


1400 
1700 
l BOO 
1900 
2000 
2100 
2200 
2300 
2400 
2500 
2400 
2700 
2B00 
2900 
3000 
3100 
3200 
3300 
3400 
3500 
3400 
3700 
3B00 
3900 
4000 


4100 

4200 


NIIMRDG=7 4300 

JNDMAX=MIN0(N,9) 4400 

GD TO (4000,5000,4000) , SAME 4500 

4000 NU( INDJC, IT)=NII( IMDIC, IT) + ) 4400 

WR I TE ( I OUT ,702 ) 4700 

1 S( SAME ) , Nll( JNDIC, IT) ,TIME,XINT 4B00 

2 ,NUMRDG, INDIC , ( E ( .1 ) , J-NUMRDG , I MDMAX ) 4900 

■ K =0 5000 

DO 4100 J = NI IMRDG , I NOMA X 5100 

K=K+ 1 5200 

U 1 NT ( NU ( INDIC , I T ) , J , I NO I C ) = X I NT ( K ) /T I ME 5300 

4100 CONTINUE 5400 

GO TO 20000 5500 

5000 TSTO( INDI C ,NOST , I T ) =TSTO ( I NO IC , MOST , I T ) + T I MF 5400 

5050 K=0 5700 

DO 5100 J =NilMR0G , I NOMA X 5B00 

K=K+1 5900 

STOINT ( INDI C , J ,NOST ) =STD!NT ( I NO IC , J ,NOST ) +X I NT ( K ) 4000 

5100 CONTINUE 4100 

WRITE! TOUT, 702)S(SAME) , MOST, TIME, X I NT , NIIMRnG , INDIC, 4200 

1 . ( E ( J ) , J=NUMRDG , I NOMAX ) 6300 

GO TO 20000 4400 

6000 TBG(NOST,IT)=TBG(MOST,IT)+TlMF 6500 

6050 K=0 6600 

DO 6100 J = Nl)MRnG , I NOMA X 6700 

K=K+1 ' 6B00 

BG I NT ( J »NOST ) =RGI NT ( j ,NDST ) +X I NT ( K ) 6900 

6100 CONTINUE 7000 

WRITE ( I OUT, 702 ) S( SAME) ,NnST,TlME,X INT.NUMRDG, INDIC, 7100 

1 <E( J),J=NUMRDG, INDMAX) 7200 

20000 IF (LAST .NE ,0 ) GO TO 21000 7300 

21000 DO 21500 K=1,NM 7400 

0021400 1=1, INDIC 7500 

DO 21200 J=1,N 7600 

GO TO ( 2 1010,2 1010, 2 101 0,21020, ? 1020, 21020,21030, 21030, 2 1030) ,J 7700 

21010 I T=i 7800 

GO TO 21040 7900 

21020 I T=2 8000 

GO TO 21040 B1O0 

21030 I T=3 8200 


22 



* 


21040 IF (TSTD(I,K,lT>.NE.O)STniNT<I,.l,K)=STniNT(l,.J,K)/TSTn(J,K,IT) 8300 

21200 CONTINUE R400 

21400 CONTINUE 8500 

21500 CONTINUE 8400 

if(nrgs.eo.o) on to 23000 0700 

00 22000 K=1,MRGS 8800 

00 21400 .1 = 1, N 8900 

J F ( B S T D ( J ) • N F » K ) 00 TO 21400 9000 

GO TO (21510,21510,21510,2152.0,21520,21 520,21530.21530,21530) ,J 9100 

21510 J T= 1 9200 

GO TO 21540 9300 

21520 JT=2 9400 

00 TO 21540 9500 

21530 I T=3 9400 

21540 BGINT ( J ,K )=R0 INT (J,K)/TBG(K,IT) 9700 

21400 CONTINUE 9800 

22000 CONTINUE 9900 

23000 INniC=INOIC-l 10000 

GO TO (25000, 25000, 25000, 23100, 23100,23100, ?3300,23300, 23300), N 1.0100 

23100 00 23200 T=1,INDJC ]0200 

23200 NU( I ,1)=MIN0(MII( I ,J ) ,NI)( 1,2) ) 10300 

on to 25000 io4oo 

23300 00 24000 I=1,IM0IC 10500 

24000 NU( I , 1 ) =M I NO ( Ml I ( ] , J ) ,NU( 1,2) ,NII( 1,3) ) 10400 

25000 CALI. MODIFY 10700 

RETURN 10R00 

701 FORMAT! ' 1 INPUT INFORMATION* /IX, 19( 1/ • DTYPF MIIMRFR TIMF ', 10900 

13 ( • INTENSITY • ) , ' NilMRDG • , 3X , ' 1N0IC * , 5X, ' FI.FMFNT* , 5X, • F| FMRN-T* . 11000 

2 5X, 'ELEMENT' ) 11 100 

702 FORMAT! 'O' ,A?,3X,J5,2X,F?.?,3T2X,FU.l ), 19, !R,RX,A?,inX,A2, JOX, A2) 11200 

1 FORMAT! F4.0 , 3FR ,0 , 1 X ,411,12,11) 11300 

7 FORMAT ( •' ' , A5 ) 11400 

, END 11500 

SUBROUTINE ENDIT (N,AK1,AK,*) 100 

DIMENSION AK1( 10) ,AK( 10) ?00 

C THIS SUBROUTINE TELLS WHFTHER OR NOT TO STOP ITERATING 300 

IC0N*0 400 

DO 34 1=1, N 500 

I F ( ABS(AK1( I )/AK( I )-l.) ,GT.5.n-3) ICnN=IC,nM+l 400 

34 CONTINUE 700 

If I ICON.EQ.O) RFTURN 1 800 

RETURN 900 

END 1000 

SUBROUTINE FLUCOR (C,FL,EL,NEL) 100 

C0MMnN/WVAL7nMEGA( 100,2 )/RVAL/RK( 100 ) /EV/ VC 1 100 , 2 ) 200 

COMMON/ XX X/X 300 

1 /INFO/ E ( 50 ) ,2 ( 50 ) ,A ( 50 ) 400 

5 /INOUT/ I OUT 500 

3 /F9/ ABCQ( 50, 50, 3), SIG, THETA 400 

4 /ARGS/ LINE(10),V0,M,L,LL,LLL,I,NCnilNT 700 

REAL*4 L0RKI2) /'K','L'/ BOO 

REAL*B W J -.8LA ( 4 ) /2*'K ALPHA ',?*'K RETA ',?*'L ALPHA •/ 900 

DIMENSION EL ( 10) ,C( 10) ,ITEST( 4) »P( 4) 1000 

DATA P/1.0,0.24,0.1,0.024,4.2,1.0/ 1100 

l)A*VO/VC ( L ,LL ) 1200 

AR«RK(L) 1300 

FL*1 .0 1400 

DO 5 J=1,M 1500 

00 5 I J=1 ,NEL . 1400 

I F( E ( J ) ,NE. EL ( I.l ) ) GO TO 5 1700 

CALL TEST( J ,L , I TEST ) 1R00 

SFK*0.0 1900 

DO 4 K=LL ,4 ,2 2000 

I F ( ITEST(K).NE.l ) GO TO 4 21.00 

GO TO (125, 125, 150, 150, 175, 175), K 2200 
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C KALPHA RADI AT.! ON 2300 

125 Kl = l 2400 

KK=l 2300 

GO TO 200 2300 

C KBETA RADIOAT JON 2700 

150 Kl=3 2R00 

KK= 1 2900 

GO TO 200 3000 

C LALPHA RAO! AT 1 DM 3100 

175 Kl»2 3200 

KK=2 3300 

200 K A Y*Z ( J ) 3400 

JFtVCUUKKI.EO.O. ) GO TO 4 3500 

IFIABCni J,J,Kl) .FQ.O.O) GO TO 4 3600 

UfUVO/VCMtKK) 3700 

DFM=0.0 3R00 

DO 3 K2=l ,M ?900 

DO 3 ,) J = l , M F I. 4000 

I F ( F(K2).MF.FL(.U) ) on TO ? 4100 

OFMsDFN + C! JJ)#ARCTHK2.J,K1 ) 4200 • 

3 CONTI MUF 4300 

XF=X/nFM , 4400 

Y=S I G/DEN 4500 

I F ( I.L .EO .2 ) AR = 4.0 4600 

FK=.5*P(K )*( J .-1 . / AR ) *OMFGA ( J , KK ) 4 ( A ( I.) /A ( 4 I ) ) * ( ( HR- 1 . ) / < IJA- 1 . ) ) 4700 

1**1. 67*1 ARCnU.J.Kl )/DEN)*( I A 1.00 fl . + XF ) /X F ) + I ALOG ( l . + Y ) / Y ) ) 4fl00 

FK=G( I J)*FK 4900 

SFK=SFK+ FK 5000 

1TP$T(K )=0 5100 

IF(NCnilNT.RO.l) WRITE! I OUT ♦ 71 21 F ( , 1 ) *K AKRI, A ( K ) » F 1 1 . ) . I.ORK { LL ) » FK * 5200 

1 FI,i),C(M) 5300 

4 CONTINUE 5400 

FL=FL+SFK 5500 

5 CONTINUE 5600 

I F ( NCOUNT . NF . 1 ) RETURN 5700 

I F ( FL.NE.l) GO TO 400 5A00 

WRITE! J OUT » 7 13) E(l.) 5900 

RETURN 6000 

400 WRITE! IOUT, 714) F(U,FL 6100 

RETURN 6200 

712 FORMAT!' « , A? , » - * , AR . • FLOURFSCPO ' , A 2 , ' - • , A ? , • . • ,5X , ' KP= * , F 1 3. 5 , 6300 

1 7X * ' C < • » A 2 ♦ ' )=• ,F13.5) 6400 

71? FORMAT!* NO ELEMENT FLOUR ESC EO • , A2 , ' . • • KF ( TOTAL ) = 1.0'. ) 6500 

714 FORMAT! • TOTAL FLOURFSCENCE CORRECTION FOR • , A2 , • = ' , F 1 3 . 5 ) 6600 

UNO 6700 


SUBROUTINE INTERP <RECU,IZ,RAB) 100 

COMMON/RA/R! 11,11) 200 

00 1 K = 1 ,10 300- 

IF! IZ.LE.K*10) GO TO 3 400 

1 CONTINUE 500 

? 1 1 = K 600 

12 = 11+1 700 

1 F ( RECU. LF . . 1 ) GO TO 7 BOO 

DO A L=2 ,10 900 

AL=L 1000 

IF(RFCU.LF.AL*.l ) GO TO fl 1100 

6 CONTINUE 1200 

7 Ji-T. 1300 

R=.01 1400 

GO TO 9 1500 

P J 1 =L 1600 

R = ( A L - 1 • ) * . 1 1700 
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9 J?*J1+1 lflOf) 

iPUi.EO.n on to in l9nn 

x*.io ?onn 

on Tn ii 2100 

10 x*.09 ??oo 

n y*io.o ?300 

RAR1=R( Il,Jl)+( f RFC1I-B ) /X ) *<« I 1 1 t.l? )-R < 1 1 ..II ) ) 2400 

RAR?=R ( I 2 »*)i J + 1 ( RFC1J-R ) / X )4( R ( J 2 , ,1? ) -R (1 2 . ,1 1 ) ) 2500 

rar=rari+( < !Z-in*(K-n )/yi*(rar 2 -rari ) 2400 

• RETURN 27 no 

FND 2 BOO 

SUBROUTINE MODIFY 100 

)NTFRFR#4 RSTD , 200 

COMMON /RACK/ RGI NT ( 9 .9 ) *RKGR ( 9 ) , RSTD( 9 ) .NRGS.RT I TLF (20.9) 300 

X /STAND/ FSTPI )0,10) ,.STniMT(9.9, lO),MFSTn| 10) , 400 

X CS( lOtlO) »STINTD(9. 10, R) ,COMPST( 9 , 10 ) 500 

X /IJNKWN/ llIMT(75,9,fl ) ,K'FII.C( 10) 400 

3 /IMFH/ F { P 5 ) /IMnilT/ I OUT 700 

4 /SCRIPT/ N.NIMfi,?) .NM. m.INDIO ROO 

5 /COR RFC/ OR! FT I 9 » 10 t B ) * TAII( R) f I SPFC ( 9 ) 900 

DFADT( X,Y)=X/( 1-X*Y) 1000 

00 4000 1 = 1 » M 1100 

RKGR I I ) =0 . 1200 

if (nrcs.ro. o) on to moo 1300 

C BACKGROUND CDRRFCT I DNS 1400 

k=rstd< i ) moo 

IF(K.FO.O) GO TO 1000 1400 

BKGR ( I) aRGIMT ( I .K) 1700 

WRlTEUOUT.i) (RTITLF(J.Ki',U=l,20).RKGR( J ) 1R00 

1 FORMAT ( • » .20A4.5X.E13.5»'' COUNTS /SECOND' 1900 

iooo jj«ispec<n 2000 

00 2000 K = 1 »NM 2.100 

NEST=NFSTD(K) 2200 

00 1100 J = 1 .NFST 2300 

I F I E( I ) .EO.ESTD(a).K) ) GO TO 1150 ?400 

1100 CONTINUE 2500 

GO TO 2000 2400 

C DRIFT CORRECTIONS 2700 

1150 00 1400 L=1 » INDIC 2R00 

DRIFT ( I »K »L ) "0 . 2900 

L 1*1+1 3000 

IF ( ABS( STDINT ('1 1 » I »K ) -STDINTI l* . I .K ) ) .GF .4.4 SORT ( STO! NT ( L . I , K ) ) ) 3100 

1DRIFT( I,K,U = (STniNT(Ll.!.K)-STDINT(U I.K) J/NIHI.,1 ) 3200 

C STANOARD INTENSITY CORRECTIONS 3300 

1500 CONTINUE 3400 

ST INTO! I t K»U=OFAOT(STOINT(L»I .K) tTAIimJ) I-RKGRU ) 3500 

1400 CONTINUE 3400 

COMPS T ( I ,K)=C5M.K) 3700 

2 FORMAT ( 31 4.9Fl?u5 1 3R00 

2000 CONTINUE 3900 

C DEAD TIME AND BACKGROUND CORRECTIONS TO UNKNOWN 4000 

00 3500 1 = 1. INDIC 4100 

N UP =NU( L 1 1 ) 4200 

00 3000 J=1,NUP 4300 

3000 UINTM,I t l-)=OEADT(UINT( J.I.U ,TAU( J,l) )-RKGR( I) 4400 

3500 CONTINUE 4500 

4000 CONTINUE 4400 

RETURN 4700 

END 4R00 


SUBROUTINE OXYGEN <OX,C,A,N,IOX,COX> inn 

dimension nxuoo) ,c( jo) ,a( 10) , 7 (io) ?no 

ct b o. 300 

nn ion i*i,w 4 no 

100 CT*CT+C( 1 ) / ( OX ( I) *A ( l ) ) *,00 

C( inx)*cnx+i 6 ,#r,T 600 

RETURN 700 

FNO 800 

SUBROUTINE TEST (KKK,LLL,ITEST) 100 

REAL*4 KALPHA ,KRETA ,KFOGF , L ALPHA ♦ L FOGF 200 

COMMON /INFO/ SPACF ( 1*50 ) ,KA(.PHA( 60 ) , KRFTA ( *50 ) ,KPOGP ( 50 ) , 300 

l L ALPHA! 50) ,1 POGF ( 50 ) 400 

DIMENSION ITFST ( 6 ) 500 

DO 100 I«] ,6 600 

100 ITFST I 1 ) “0 TOO 

J F( KALPHA ( KKK ) . IT. K EDGE (I.LL ) ) ITFST ( ] ) «1 000 

IF(KALPHA(KKK) .I.T.LFOGF(LLL) ) ITFST(2)*1 . 900 

I F( KRFTA ( KKK ) . I. T.KEDGE ( LLL ) ) ITPST(3)»l 1000 

IF (KRFTA (KKK).LT.LFOGF(LLL) ) TTFST(4) = ). 1100 

I F( L ALPHA ( KKK ) .LT. KFOGF ( LLL ) ) ITEST(6)«l 1 £00 

I F (LALPHA ( KKK ) .LT.LEDGE(LLL) ) ITFST(6)*1 1300 

RETURN 1400 

END 180C 

BLOCK DATA 100 

COMMON /RA/R( 11,11) 200 

COMMON/SS/S < ? ) 300 

OATA R/l.00,.934,.856,.786,.735* .093 , .662 , , 635 , .411 , .59? , .578 , 400 

Al.00,.944,.R73,.R0P,.760,.71R,.68R,.663, . 630 , . 6 1 ? , . 606 , 1 .00, .953, 500 

B. 88R,.fl*2fl,.7fl2,.741,.713».687,.665». 639 ,.634,1 .00, .961, .903, .847, 600 

C. R04, .764 , • 737, .71 3 » .691 , .665 , . 661 , 1. 00, . 968 , .917, » 867 » • 827, • 789, 7Q0 

0. 764,. 740.. 71R.. 695,. 691 ♦ 1 .00 . .975 , . 933 , .88R,.B51 . .fl 1 7 , .793 , .770 , ROD 

E. 750, . 730 , .725, 1 .00, .981 • .948 , .91 1 , .878, .847 » • 825, • 805, .785, • 767 ♦ 900 

F. 763,1. 00,. 9flfl,. 963,. 935,. 907,. RRl , .862 , .844, ,826 , .fill , .806 , 1 .00, 1000 

G. 993,. 977, .959,. 93R,. 919,. 904,. 889,. 874,. 862,. 858. 1.00,. 997,. 990, 1100 

H. 981, .970, .9 59, ,9 50, .941, .932, .924, .921, 11*1. 00/ 1200 

DATA S/'U* ,’S» ,*8'/ 1300 

ENO 1400 
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APPENDIX II 

Glossary of Terms Used in the Microprobe Correction Program 


A(I), 1- 1 - M 
AK(I), 1^I<N^M 


AK1 (I), 1<I<M 


Atomic weight of the element, E(I). 

Ratio of the measured intensity readings for E(I) 
in the unknown to the average corrected intensity 
reading of E(I) in standard K. (Equation 2). 

Value of k A (calculated) at the ith iteration; 
AK1(I) * AAA* UC(I)/SC(I). 


AAA 


AB 


ABN 


ABDS(I, K), 1<I<N 

1< K<NM 

ABS 

AT 


ATN 

ATDS(K) , 1< K<NM 


ANC 


ABC0(I, J, K) 

I ^ M 
1 £ J £ N 
1 ^ K< 3 


Product of ZAF. This is calculated separately 
for each element E(I). 

Dummy argument for absorption correction in 
subroutine ABSCOB. 

Absorption correction for the unknown. This is 
computed for each element E(I). 

Absorption correction for E(I) in standard K. 


Combined absorption correction 
ABS = ABN/ABDS(I, K). 

Dummy argument for the atomic number correction 
in subroutine ATNCOR. 

Atomic number correction for the unknown. 

Atomic number correction for standard K. 

Combined atomic number correction 
ANC « ATN/ATDS(K) 

Matrix of mass absorption coefficients with E(I) 
the absorber, E(J) the emitter, at the K a -line 
(K = 1), Kp -Une (K = 2) or L a -Une (K = 3). 
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BSTD(I), 1 ^ I - N If BSTD(I) = K, 1 £ 9, then background standard 

K is used for E(I), If BSTD (I) * no bac kground 
correction is necessary. 

BKGI<(I), 1 ~ I ~ N Average background intensity for E(I), 

BGINT(I, K) 1 S I & N Average background intensity readings for each 

IS KS NBGS element E(I) using background standard K in units 

of counts per second. 

B TITLE (J, K), J-1,2 ... 20 Identification card for background standard K. 

1 < K < NBGS 

0(1), 1SI<NEU 


CK(I) 

CKL(I) 

CL1(I) 

CL2(I) 

CLM(I) IS I <M 
CMX(I) 

CM2(I) 

CM3 (I) 

CMN(I) 

CS(I,K), 1SISNESTD (K) Composition (wt%) of each element (ESTD(I,K)) 

1 S K S NM present in standard K. 

C0XY Wt% oxygen present in the unknown which is given 

as input, C0X is dummy argument for C0XY', in 
subroutine OXYGEN* 

C0MPST(I, K),1<I< N Wt% of E(I) present in standard K. 

1< K<NM 

CHI X) - used in computing absorption correction 

(Equation 8, 9). 

DRIFT(I,K,L),1<I<N 
1<K<NM 
lS L<INDIC-1 


Total Drift correction to intensity readings for 
E(I) in standard K* L refers to INDIC-1, the 
number of times the standards have been 
measured. 


Composition (wt%) of E(I) present in the unknown. 
C is used as a dummy argument denoting the 
composition of a sample in subroutines FLUCOR, 
ATNCOR, ABSCOR, OXYGEN 


Constants for computing mass absorption coeffi- 
cients in subroutine ABC0 where p/p sC K n . 
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EL 


DR 

DATA 


m 


ESTD(I,K) 

1< I < NESTD(K) 
1* K^NM 

FL 


FLN 


FLDS(I,K), 

1< I <N 
l£ K*NM 

FLU 


IN 

I0UT 

10X1 


Dummy argument used to specify chemical symbol 
of each element in a sample. EL is used in sub- 
routines ABSCOR, FLUCOR and ATNCOR. 

Drift correction for unknown calculated for each 
element E(I). 

NAMELIST name which inputs the following param- 
eters: THETA, TAU, STD, N, NEU, M0DE, VO, 
LINE, ISPEC , NPR0B, NSTC0M, NM, NBGS, BSTD. 

Array of chemical symbols of all elements used in 
the problem; this table is ordered as described in 
Appendix IV (E File). 

Array of chemical symbols for the elements E(I) 
present in the standard K. 


Dummy argument for fluorescence correction in 
subroutine FLUCOR. 

Fluorescence correction for the unknown. If nec- 
essary, this is computed for each element E(I) in 
the unknown. 

Fluorescence correction (if any) for E(I) in 
standard K. 


Combined fluorescence correction 
FLU = FLN/FLDS(I,K) 

FLU = 1.0, if no correction. 

Specifies which computer input unit is used. 

Specifies which computer output unit is used. 

Subscript used to denote whether oxygen is present 
in the problem. 

If 10 X 1 = 0, no oxygen 
If 10X1 > 0, oxygen is present 
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ISPEC(I), 

1= 1,N 

INDIC 

J1(I), 1 = 1,M 
J 

KEDGE (I), I = 1,M 
KALPHA(I), I = 1,M 
KBETA(I), I = 1,M 
LALPHA(I),I = 1,M 
M 

N 

L1(I), I = 1,M 
L2(I), 1= 1,M 
L3(I), I - 1,M 

M1(I), M2(I), M3 (I), 
M4(I), M5(I), I = 1,M 

NK(I), I = 1,M 
NKL(I), I = 1,M 


Denotes whoch spectrometer was used to measure 
E(I). 

The number of times the standards as a group 
were measured (see Appendix IV, Data set D). 

Mean Ionization Potential - Input to MAIN. 

Mean Ionization Potential - as used in subroutine 
ATNCOR. 

K-edge for E(I) 

Wavelength of K a -line for E(I). 

Wavelength of Kg --line for E(I). 

Wavelength of L a -line for E(I). 

The number of different elements used in any 
problem, 1< M^ 50. This includes elements in 
standards, and unknowns. 


Number of elements in the unknown measured by 
the probe. 

Wavelengths associated with the 3 L-edges, 

L r , L ri , L m ; used in subroutine ABC0 to obtain 
mass absorption coefficients. 

Wavelengths associated with the 5 M-edges , used 
in subroutine ABC0 to obtain mass absorption 
coefficients. 

Exponent of K for each element (I) in the expression 
for mass absorption coefficient (n/p = C \ n ) in the 
wavelength region from 0.7A to the K-edge of ele- 
ment I. 

Exponent of k for each element (I) in the expression 
for mass absorption coefficient (p/p = Ck n ) in the 
wavelength region from the K-edge to the L ni -edge 
of element I. 
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NI(I),I ~ 1, M 

NEU 

NM 

NU 

NC 

NBGS 

NESTD(K) 

K = 1, NM 

NEST 

NSTC0M 

0XY(I), I = 1,M 

0MEGA(I,K), , 

1 = 1, M 
K= 1, 2 

RB 

BBS(K), 1< K< NM 
RBU 

RK(I), 1 < I < M 

R(I,J), 1=1,11 
J = 1, 11 


Wavelength associated with the first N-edge; used 
in subroutine ABC0 to obtain mass absorption 
coefficients. 

Total number of elements in the unknown. 

Number of standards used in a problem. 

Number of readings taken on an unknown sample. 

NU is obtained by the computer internally. 

The number of the unknown sample being analyzed. 

NC is obtained by the computer internally; 1 - NC - NU. 

Specifies the number of background standards used 
in the analysis. 

Number of elements in standard K 


NEST = NESTD(K) 

See Explanation of input data set A (Appendix IV). 

0XY(I) = N/M if the oxide is of the form E(I) N 0 M . 

Fluorescence yield ft (I, 1) for K radiation, Fluo- 
rescene yield ft (I, 2) for L radiation. 


Dummy argument for backscatter coefficient in 
ATNCOR. 

Backscatter coefficient for standard K. 

Backscatter coefficient for unknown. 

Absorption jump ratio, r, for K lines in each 
element E (I). 

Tabulated values of R. I indicates the atomic 
number and J indicates the reciprocal of U = 
VC/VO. 


*! 


RECU 


Reciprocal of U * VC/VO. 


STDINT (J, I, K), 
J = 1, INDIC 

X ^ I * N 
1<K<NM 


Intensity of each element E(I) in standard K 
measured in counts per second. J gives the num 
ber of times the standards as a whole have been 
measured. 


STINTD(I, K, L), 
1^ I* N 
1< K < NM 
1< L< INDIC-1 


Intensity readings in counts per second for each 
element E(I) in standard K after correction for 
dead time and background. 


STITLE(K, J) Identification card for standard K. See Appendix 

1 < K < NM IV, Input file C. 

J = 1, 20 


SIGMA er (Equation 11). 

SB Dummy argument for stopping power in ATNCOR. 


SBS(K), 1 < K < NM Calculated stopping power for each standard K. 

SBU Calculated stopping power for each unknown. 

STD(I, J) See Appendix IV, input data set A. 

1< I< M 

l-<- NSTC0M 

STINT The fully corrected intensity for each element in 

a standard. 

SCON(I, K) Weight % of element E(I) present in standard K. 

1 < I < M 
1< K< NM 

TACJ(J), J =1,3 Dead time for spectrometers 1, 2, and 3 

respectively. 


TBG(K, I), Counting time on background standard K. 

K = 1, NBGS I indicates which group of elements is 

1 < I < 3 being measured. 
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TSTD(I, K, J) Counting time on standard K. I gives the number of 

1*1* INDIC times all of the standards have been measured. J 

K = 1, NM indicates which group of elements is being measured. 

J = 1, 3 


THETA Take off angle of the electron probe. 

UINT(NC, J, L) Intensity reading for an element E(J) in counts per 

1 < NC - NU(L) second in the unknown. NC indicates the number of 

1 ~ <3 - N the sample and L indicates the number of times the 

1- L- INDIC -1 standards as a whole have been measured. 

UTITLE (J), Identification card for the unknown. 

J *= 1, 20 

UC(J), 1< J< N Calculated wt. % of element E(J) present in the 

unknown. 

VO Operating voltage. 

VC (I, LL) Excitation energy of each element E(I) for either 

1 £ I £ M the K-line (LL = 1) or the L-line (LL = 2). 

LL = 1, 2 

X Dummy variable for CHI in subroutine F LUC OR 

and ABSCOR 

Z(I), 1< I< M Atomic number of E(I). 
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APPENDIX III 

Detailed Descriptions of Subroutines 


A. MAIN - This subroutine inputs the data files (A,B,C) and calls DATPTS to 
input data set (D). For those elements whose compositions are specified in the 
unknown, it is assumed that these compositions will not vary from unknown 
sample to unknown sample. MAIN calls for the calculation of mass absorption 
coefficients, and for each standard calls for the computation of atomic number, 
absorption and fluorescence corrections . To each element in the unknown MAIN 
assigns a standard and then calculates in turn an estimate for the composition 
of the unknown, the corrections associated with the estimated composition and 
then re-estimates the composition. This is done for each set of intensity read- 
ings on the unknown. If MODE = 2, the composition of the sample is known, and 
the intensity ratios of unknown to standard are calculated for each element, 
Finally, the MAIN routine outputs a detailed account of the computations to assist 
the experimenter in his analysis. Figure 1 presents a flow chart of MAIN. 

B. DATPTS - DATPTS, as developed for the GSFC microprobe, reads and 
processes data in the following fashion: Each card, which DATPTS inputs, con- 
tains a time and 3 intensity counts, one for each spectrometer, as well as other 
information outlined for the Intensity Data File (Appendix IV). 

If the readings are taken on the sample, the intensities are divided by the 
time to obtain the number counts per second for each measured element. The 
total number of data points taken on the sample are then computed. If the read- 
ings are taken on the standards or backgrounds, the average intensities and 
average counts per second are obtained. 

When all the readings for the problem have been input, DATPTS calls sub- 
routine MODIFY which computes corrections to the raw data. The flow chart of 
DATPTS is shown in Figure 2. 

C. MODIFY - To each measured element MODIFY assigns a background correc- 
tion according to the values of NBGS and the BSTD array (Appendix IV). It com- 
putes the drift correction for each element standard, storing the corrections in 
the DRIFT matrix. The count rates for each standard are corrected for dead- 
time and background and stored in the STINTD matrix. Then the unknown sample 
count rates are corrected for dead-time and background and restored in the 
UINT matrix. Figure 3 illustrates the flow of MODIFY. 
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D. ATNCOR - for each sample this subroutine computes the atomic number 
correction according to Equations 6-7. The (dummy) arguments of the sub- 
routine are: 

RB - the backs cutter coefficient 

SB - the stopping power 

NEL **• the number of elements in the sample 

AT - the atomic number correction 

EL - the element array for the sample. 

Figure 4 gives a flow chart of ATNCOR. 

E. ABSCOR - computes the absorption correction for each element in each 
sample (Equations 8-11). NEL, C, EL, AB are the arguments of the subroutine, 
ABSCOR. The first 3 have the same meaning as in ATNCOR and AB is the ab- 
sorption correction for the element E(I). The MAIN program selects E(I) and 
transmits the parameter, I, in the ARGS COMMON block. Figure 5 illustrates 
the flow of ABSCOR. 


F* FLUCOR - computes the fluorescence correction (if any) for each element 
E(L) in each sample (Equations 12-14). The arguments EL, NEL, C are the same 
as in ATNCOR; FL is the fluorescence correction and L is transmitted in the 
ARGS COMMON block,. Figure 6 is a flow chart of FLUCOR. 


G. TEST — is a subroutine called by FLUCOR to determine if an element, E(J), 
in a sample has fluoresced the element, E(L), measured for that sample. When- 
ever fluorescence occurs, the ITEST argument is set to 1; otherwise it remains 
0: 


If K a E(J) < K e E(L), ITEST (1) = 1 

K a E(J) < L E(L), ITEST (2) = 1 

K E(J) < K e E(L), ITEST (3) = 1 

K E(J) < K edge E(L), ITEST (4) = 1 

lf a E(J) < K edge E(L), ITEST (5) = 1 

L a E(J) < L edge E(L), ITEST (6) = 1 

J and L are also transmitted as arguments. 


H. INTERP - linearly interpolates from the tabulated values for R given by 
Duncumb and Reed. 4 This subroutine is used by ATNCOR to compute the value 
of the backscatter coefficient. The arguments are IZ, the atomic number of E(I); 
RECU, (l/U); and RAB, the value of R. The values of R 4 are contained in the 
block data at the end of MAIN. 
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I. ENDIT - Tests for convergence for each measured element, E(I) in the un- 
known sample. If sufficient convergence is achieved, Equation 18, then ENDIT 
directs MAIN to end the iterative procedure and output the final composition 
of the unknown. If not, ENDIT returns to the iterative portion of MAIN. 
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APPENDIX XV 
Description of Input Files 


A-DATA FILE ~ The Fortran IV namelist feature directs the input of DATA 
File. This option enables data to be identified by name; for example, 1 = 1 will 
cause a 1 to be placed in the location labeled I; E(3) = 'Fe' causes the third item 
in the E-array to be identified as iron. The parameters initialized in this 
fashion are described below. The statement in the MAIN program 


100 READ (IN, DATA, END = 10000) 


accomplishes the initialization. Figure 8 shows an example of the sample 
NAMELIST cards for one of the problems presented in Appendix VI. 

M is the total number of elements present in both the unknown and standards 
NEU is the total number of elements in the unknown 
N is the number of elements measured by the probe 
NM is the number of standards used in an analysis 
NSTCOM is the number of ways in which the standards were combined to 
determine the composition or intensity ratios for the unknown. This allows 
different standards to be used for calculation of the same element. For example, 
if element A appears in both standard #1 and standard #2, the experimenter may 
wish to use standard #1 first for A to obtain the amount of A present and then 
standard #2 for the same purpose. Up to 10 such combinations may be made 
for a problem. 

STD(I,J), I = 1, N ; J = 1, NSTCOM specifies which standard is to be used for 
E(I) in standard combination J; for example if STD (4,2) = 3, use standard #3 
for E(4) in the second analysis of the data. The standard is designated by number 
which corresponds to its order in the Composition File. 

NBGS specifies the number of background standards used in the analysis : 

0 < NBGS < 9. This program assumes one background standard per element. 

BSTD(I), I = 1, N - If NBGS > 0 there is at least 1 element measured for 
which the background is non-zero; hence, the BSTD array must be initialized for 

1 = 1,2, . . ., N, in the following manner. If the background for E(I) is zero, set 
BSTD(I) = 0; if the background is not 0, set BSTD(I) = K where K is the number 
of the background standard used for E(I). 
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MODE - If MODE « 1, intensity readings from the probe (Intensity Data 
File) are input and the weight % present in the unknown is determined for the 
measured elements. If MODE « 2 the composition of the unknown is known and 
the output consists of the intensity ratios of unknown intensity to standard in- 
tensity for each element in the unknown. 

UNE(I), I = 1, M specifies whether the element, E(I), was measured with the 
K or L x-ray line. 

ISPEC(I), 1*1, N specifies on which spectrometer (first, second, or third, 
etc.) the element, E(I), was measured. 

VO gives the operating voltage for the probe (E 0 ) in kilovolts. 

TAU(I), I = 1, 2, 3 specifies the dead -times for spectrometers, 1,2,3 
respectively. 

Note: If the probe has more than 3 spectrometers, the dimension on TAU 
should be increased accordingly. 

THETA is the take-off angle given in degrees. 

B- Element File 

For elements 3-92, file cards are provided for each element. Each of these 
sets of 3 cards contains the element’s chemical symbol, atomic number (Z), 
atomic weight (A), excitation energies for the K (elements Z = 5-39) and L (ele- 
ments Z - 12-92) lines in kilovolts (Bearden), 18 mean ionization potential (J*) 
in electron volts (Duncumb and Reed), 4 the K absorption jump ratio (r A ) (Colby), 19 
the cation to anion ratio if oxide is formed (OX), the K shell fluorescence yield 
(W(B)) and the mean L shell fluorescence yield (W(B)) (Fink, et al.). 20 the K a 
and wavelengths for elements Z = 3-39 and the L a wavelengths for elements 
Z * 20-92. For elements less than atomic number 20, the L a wavelength is set 
to 50A. Also given in the E file are the K, Lj» L n , L iri , Mj , M n , Mjjj , M IV ,, 
M v ,Ni wavelength edges (Bearden) 18 in A units and the CK, NK, CKL, NKL, 

CL1, CL2, CLM, CM1, CM2, CM3, CMN parameters used by Henrich 6 to calcu- 
late the mass absorption coefficients. 

To enable the user to calculate approximate mass absorption coefficients for 
L Hnes of elements of atomic number 24-29, values of CL1 and CL2 were de- 
termined by linear extrapolation from the values given by Heinrich 6 for elements 
Z = 30 and above. The complete Element File compiled to date is given in the 
following Table. 
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TABLE I 
Element File 


FL NO. WGT VCf K >- VC<LI- 

-NK- KFDG6 CKL- NKL- 

-CMl- -M? CM? M3- 

•LI* 3. 007.0 0.054R 00.000 
?.RB ?26.5 

' Rfi • 4. 009.0 00. in 00.000 
?.R6 111. 00 

• R* 5. 010. P 00.000 00.000 

?.R5 

• C' 6. 01?. 0 00.2R4 OO.O0O 

?.P* 43.f.R0 

• N* 7. 014.0 00.400 00.000 

?.R3 30.990 

• n* fl. 016.0 00. *53? 00.000 

?.P? ??.??0 

• F • 9. 0)0.0 00. AO? 00.000 

?.fll 

•MFHO. 0?0.? 00.867 00.000 
?.P0 14.300 

•NA'l). 0?3.0 01.070 00.000 
? .70 11.480 O.A? ? .73 

•MG*)?. 0?4 .3 01,300 00.049 
?.79 OQ. 41? 000,89 ?.73 

•ALM3. 0?7 .0 01.560 00,073 
?.7P 07.951 001 . ) 8 ? .73 

• S f '14. 028.1 01.840 00.100 
? . 7 7 -06.738 001,54 ?,7? 

• P • 1 5. 031.0 0? .140 00.13? 

? .76 05.7R4 00?, 43 ?.73 

• S * 16 . 03?. 1 0? .470 00.000 

? ,76 OA.Olfl 00?. 43 ?* 7? 

• GUI 7. 035.5 0?.fl?0 00.000 

? .76 04.377 00?. 9ft ? .73 

• AR • IR • 039.9 03. ?00 00.000 
? • 7 5 03.P71 003. A? ?,7? 

• K' 19. 039.1 03.610 00.?95 

? .75 03,437 004.31 ?.73 

•CA*?0. 040.1 04.040 00.349 
?,74 03.070 005.10 ?.7? 

•SC*?1. 045.1 04.490 00.000 
3.74 0?. 7620005 .99 ?.73 


-J~ -RK— -OXY-- OMGAK KALPHA KRFTA- (.ALPHA — CK— 

—11 CL1 L2-- -CL?— — L3 CLM— —Ml — 

■ CM3— — *M4— — — M5— — f.MN— — Nl— — OMRAI., 

0*000 00.00 1.0000 .0000 00.000 00.000 50.000 000.135 

,00375 

0000 00.00 1.0000 ,0000 00.000 00.000 50.000 000.360 

.00500 

0000 15.73 0.6667 ,0000 A7.A00 00.000 50.000 000,740 

,00626 

0146 15.14 1.0000 .0009 44,700 00.000 50.000 001,360 


.00750 

Ol?5 14.59 1.0000 0,00? 31.AOO 00.000 50.000 00?.?10 

. OOP76 

Oi?7 14,08 1.0000 0.004 ?3.700 00.000 60.000 003. POO 


,01000 

01 ? 3 13. AO 1.0000 0.007 18.300 00.000 60.000 004,000 


' . 01 1 ? 6 

01 ?3 13.16 1.0000 0.015 14.600 14.600 50,000 006,770 

.01 ?50 

01? A 1 ?. 74 2.0000 0.0?0 11.910 11.6)7 00.000 009.060 


,01.375 

0133 l?. 35 1.0000 0.030 09.889 09.668 60.000 01.) .760 

?60. 7 

.0)500 

0)^? 11.99 0.6667 0.04? 08.337 07.981 60.000 014.870 

170.4 

.01 A ? 5 

0)64 11.64 0.5000 0.06A 07.1?6 06.768 60.000 018. 500 

1?3.0 

.01750 

01 AA 11.31 0.4000 0.075 06.156 06.804 60.000 0??.«00 

96.0 

.01876 

0180 11.00 1.0000 0.096 05.373 05.03? 50.000 027.000 

, 0?000 

0194 10.71 1.0000 0,1.17 04.7?8 04.403 60.000 031,700 


,o?l?5 

0209 10.43 1.0000 0.138 04.193 03.886 50.000 036.900 


,0??50 

02?4 10.17 ?. 0000 0.159 03.74? 03.6-54 50.000 042.600 

42,1 

.02375 

0239 09.92 1.0000 0.180 03.359 03.090 36.330 068.400 

35.5 

.02500 

0255 09. AR 1.0000 0.?0? 03.03? 0?.7R0 31.350 065.100 

. 0?A?5 
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TABLE I-(continued) 


PI. Mfl, WGT.- VCIKl- VC(U- -J-- -RK— -flXY — flMOAK KALPHA KRFTA- LAI, PH A — CK — 
-NK- KPnPF- — CKl - MKL- —LI — "CL!— — L?— -CL?-- — L3— -CLM- — M l — 
-M2— CM?- —M3— f,M3- — M4— —— M*5— — f.MN— —Ml— DMCA I, 

•TH22. 047.9 04.970 00.450 0?70 09.45 1.0000 0.2?4 02.748 02.514 ?7,4?0 042.100 
?,73 02.497 007.00 ?,7? ?7.3 

,0?75O 

' V* ?3. 051.0 05.470 00.000 02P6 09. ?4 1.0000 0.?46 02.504 02.284 24.2*0 069.800 
? .73 0?.?69 008.0? 2.73 

,0?R75 

•CR'?4. 05?. 0 05.990 00.598 0301 09.03 1 .0000 0.?6R 0?.?91 02.0B5 ?1*640 078,000 
?.73 02.070 009. JR ?.7? 14.7 7.R0 17.9 4.8 ?0.7 

,030on 

• MM • ? 5 . 054.9 04.540 00.000 03)4 0R.R3 1.0000 0.?90 0?.l03 01.9)0 19,400 ORA. 700 

? ,7? 01.R94 0)0.45 ?,73 R.R5 5.9 19.1 

,031 ? 5 

•FF'?4. 055.9 07.1)0 00,707 033? OR. 44 1.0000 0.3?0 01.937 01.757 17.A00 095. ROO 
?.7? 01.743 011.75 2.7? 17.2 * 10.0 17.? 7.0 17.525 

, 03? 50 

•C0‘?7. 050,9 07.710 00.779 034? 08.46 1.0000 0.?50 01.790 0).4?l 14.000 10A.400 
? .7! 01.400 013. ?5 ?.73 15.41 1 t.?5 15*41 R.0 15,9? 

.03375 

*MI*?R. 050,9 00.330 00,054 034? 00. ?9 1 ,0000 0.3R0 0J..A59 01.500 14.400 115.900 
?, 71 01.400 0)4.00 ?,7? 14.2 12.4 14.? 9,0 14.5? 

.03500 

*f,H'?9. 04? . 5 00.9P0 00.03? 0377 OR. 1 ? 1 .0000 0.410 01.54? 0) .30? 13,400 124,800 
? .71 01.3RO 014.45 ? .73 ll.?70 014. OA 13.0) 10.1 13. ?9 

, 0 3 A ? 5 

« 7M*30. 0A5.4 09.AA0 01.030 039? 07. 9A 1 .0000 0.440 01,43/, 0) ,395 12.300 13R.OOA 
?.7'0 01 . ?R3 01 R . ?5 ?.??- 10. ??001 5,40 ll.P-70011.?0 1?.1 3 

,03750 

T,A*?1 . 049.7 10.400 01.1?0 0407 07. R). 1.0000 0.4R/I 01,341 01.?0R 11.3)3 149,5 00 
? , 70 01.194 0?0.?0 2.7309.520 017. ?5 10.R30 13.40 11.100 

,03575 

•OP' 3?. 07?. 4 11.100 01,220 04?? 07.44 1.0000 0.5?0 0) ,?55 01 .1 ?R 10.45A 1 A?,?00 
? . 70 01.117 0?? . ) 5 ?,7? 0R.7?0 01R.95 09,940 0)3. AO 10.190 

.04000 

' A S • ? 3 . 074.5 I). 900 O1.??0 04?7 >07.51 1.0000 0.5A0 01.177 01.05A 09.A71 175,400 
? , A9 01.045 0?4.?5 ? . 73 OR. 107 ?0.750 09.1?4 014.90 09.390 004. ?R 

.041 ?5 

»SF»34, 079.0 1? .400 01.430 045) 07. 3R 1.0000 O.AOO 01.10A 00.990 03.990 189.400 
?.49 0O.9R0 0?4 . 40 ?.7? 07.504 ??.550 0R.4JA 1A.?00 0R.A70 004. A9 

,04? 50 

• RR ' 35. 079.9 1,3.500 01.550 0444 07. ?4 1.0000 0.A4O 01 .041 00.931 OR. 375 ?05.000 

? ,49 00.920 02R.R0 ?.73 04.970 ?4.400 007. RO 17.700 OOR. 00 005,13 

,04375 

•KR*34. OR3.7 14.300 01.4R0 04R1 07.1? 1.0000 0.4AA 00.931 00.R7A 00.000 319.300 
?.4R O0.P44 031.25 ?.7? 6.4400 24.700 007.?) 019.15 007.4? 005. AO 

.04500 

•Rfi‘37. OR5.5 1 5. ?00 01.810 04«5 00.00 1.0000 0.A9? 00.9?7 00.827 07. 31 A ?35.300 
? , Afl 00.R14 033.90 ? *73 05.99R 02R.R0 0A.A43 030.R0 004. R9 006,10 

.04 A? 5 

•5R'3R. 0R7.6 14.110 01.940 0510 00.00 1.0000 0.71R 00.R76 00.7R1 0A.3A3 351.300 
?.4R 00.770 034.50 3.73 05.583 031.30 0A.17? 033.40 OA.3R7 00A.A3 

.04750 

• Y'39. 0RR.9 17.040 0?.0R0 05?4 00. 001. 0000 0.744 00.R3O 00.73R 04.449 3AR.100 

? .67 00.7?fi 039.30 2.73 05.232 033.50 05.755 0?4.10 05.9AI 007. ).R 

. 04R75 

• ZR »40. 091 .2 00.000 02.220 053B 00.00 0.5000 0.770 00.7R7 00.700 06.070 

042.30 2,73 04.867 036.15 05. 37R 0?5.9,5 05.5.3R 007. 7fl 

.05000 

• NR • 41 . 092.9 00.000 02.370 0553 00.00 1.0000 0.796 00.746 0<5.7?4 

045.50 2.73 04.5R1 038.90 05.026 037.90 05.333 008.39 

.05400 
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TABLE I-(continued) 


FL NO . 
-NK- 
-CM1 

wot.- vo ini- vcai- 

, KFDPF CKL- NKL- 

- -M7— - CM?- — M3-" 

» wn • 4? . 

096.0 

00.000 07.520 
048.50 2.72 

•TC'43. 

099.0 

00.000 02.680 
052.10 2.72 

• RIM 44. 

101,7 

00.000 07.840 
055.60 ?,7? 

IRH'45, 

107.9 

00.000 03.000 
059.30 7.7? 

•PD'46, 

106.4 

00.000 03.170 
063.J0 ?.7? 

•60*47. 

) 07 . 9 

00.000 03.350 
067.30 ?.7i 

•CO *40. 

1 17.4 

00.000 03.540 
071.00 ?.71 

• I M *40, 

1 14,0 

00.000 03.730 
075.50 7.71 

•5M150, 

110.7 

00,000 03.930 
070,60 7.71 

*50*61 » 

1 7 1 . F 

♦ 

00.000 04.J30 
004.70 ?,70 

* TF *57. 

1?7.6 

00.000 04.340 
000.50 7.70 

* 1*6?. 

1 ?6.o 

00.000 04.560 
094.00 ?,70 

1 XF *54. 

131.3 

no. non 04.700 

090.30 7.6° 

•0.5*65, 

13?,9 

00.000 05.01,0 


•H*?"— -CLM* 


•Ml — 


>67 00.00 1.0000 0.072 00.710 05,406 

04R.50 7.7? 04. 298 041.45 04.710 079.73 04.9)3 004.05 

.06? 

>01 00.00 1.0000 0.040 00.000 

05?. 10 ?.7? 04.060 044.50 04.436 031.95 04,63? 009.74 

.060 

>95 00.00 1.0000 0.R74 04.046 

055.60 ?,7? 003.83 047.50 04.180 034.10 04.369 010,47 

,074 

>09 00,00 1.0000 0,900 * 04.597 

059.30 ?,7? 03.626 050. 70' 03.94? 036.40 04,130 011.24 

.000 

»?3 00.00 1.0000 0 ,9?6 04,360 

063.10 ?. 7? 03.4?0 063.96 03.724 0?0.70 03.900 012.03 

ii .006 

.37 00.00 2.0000 0.96? 04,164 

067.30 ?,7i 03.?54 057.50 03.514 041*30 03.690 ni?.07 

.09? 

i5l 00.00 1.0000 0.970 03.956 

071.00 ?. 71 03.005 060.70 03.3?6 043.60 03,504 013,76 

.n«o 

,65 00,00 1.0000 1.004 03.77? 

075.50 ? ,7J 0?.9?6 064.50 03.147 046.30 03.374 014,60 

, 104 

79 00.00 1,0000 1.030 03,600 

079.60 ? . 7 1 0?,777 060.00 0?. 90? 040,00 03,166 015.66 

.110 

,9? 00.00 0.6667 1.056 03,439 

004. ?0 ? . 70 0? .639 071.96 03.030 051.70 3.0000 0)6.66 

.116 

06 00,00 i.oono l.on? 03.700 

000.60 ?.70 0?. 5) 1 076.66 0?.6R7 054.30 0?.066 017.7? 

.17? 

?o oo.oo i.onon l.mo o?.l4o 

094.00 ? , 70 0? . 309 000.35 07.553 057.70 07.719 0)0,0? 

. 1 ?0 

34 00.00 1.0000 1.134 00.000 

090.30 ?.6<> 0?.?74 004.00 0?.4?9 060.30 0?.6«? 019.97 

.134 

'47 00.00 1.0000 J.160 0? . 00? 

103.40 ?.69 0? . 1 67 000.40 07.314 063.40 07,474 0?1 .16 10.20? 

0 IP.? .140 

• A A ' 56 . 137.4 00.000 05.P50 0761 00.00 1,0000 1.106 Op.775 

109.00 ?.66 0?.060 093.15 0?.?04 066.90 0?. 363 0??.41 09.067 

019.3 10.045 16.0 .146 

•LA* 67. 130.9 00.000 05.400 0775 00.00 1.0000 l.?1? 0P.666 

114.40 ? ,60 01.973 097.00 07.10,3 070.70 07,750 n?3.70 09,04? 

070.4 10.371 16.9 .15? 

• C F ' 5 0 . 140.1 00.000 05.7? 0700 00.00 0,6667 1.730 07.661 

170.00 ?.60 01.009 .102.50 07.01 1 073.60 07.164 0?5.05 00,614 

071.6 09.692 17.9 10.400 15.5 .160 

• PR • 59. 140.9 00.000 05.960 0802 00.00 1.0000 1.264 07.463 

125.50 2.68 01.811 107.30 01.924 77.000 02.077 076.47 nn.iOO 

0??.0 09.P50 10.9 09.957 16.? .164 

•MO '60. 144.3 00.000 06.210 0015 00.00 l .0000 1.290 0?,370 

132.00 2.67 01.735 112.00 01.043 01.000 01.995 077.93 07.031 

074.1 00.773 70.0 09.499 17.2 .170 

* PM '61. 145.0 00.000 06.460 0879 00.00 1.0000 1.316 07.703 

130.50 2.67 01,665 110.50 01.76? 005.00 01.910 079.44 07.513 

075.4 00.376 71.0 09.115 10.7 .176 
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TABLE I-(continued) 


FL NO. WGT.- VC(K)- VC(U- -J RK— -DXY— OMGAK KALPHA KRFTA- I.AI.PHA — CK- 

-NK- KFOGR CKL- NKL -LI CLI IP CL? L3 CLM Ml — 

-CM l M?- — CM? M3 — CM3- — M4 — -MS- — CMM — — Ml — DMGAI. 


5^*62. 150.4 00.000 04.720 0843 00.00 1.0000 1.342 

143.50 2.67 01,599 122.50 01,703 088.00 
026.7 OR. 023 22.2 OR. 705 19.1 

152.0 00.000 06.980 0856 00,00 1.0000 1.36R 

150.00 2.66 01.536 128.00 01.626 092.00 
02P.1 07.642 23.3 OR.3’1 20.1 10.604 10.893 13.05 
0*64. 156,9 00.000 07.240 0870 00.00 1.0000 1.394 

01.477 134.00 01.561 096.30 
21.2 10.111 10.40R 14.30 
0RP3 00.00 1.0000. 1,420 
01.421 140.00 01.501 100.30 

14.80 

00.000 07.790 0R97 


511*63. 


00.000 

157.00 2.66 
029.6 07.322 24.5 07.997 
R • 65. 15R.9 00.000 07.520 
163.50 2.65 

.1 07.011 25. fl 07.6P5 22.3 09.69? 09.957 
’•66. 162.5 


031 


035.9 06.179 
'69. ]. 6 P , 9 00,000 


00.00 1.0000 1.446 

170.50 2.65 01.365 146.00 01.438 104.50 

032.6 06.715 27.0 07.383 23.4 09.2*58 09.5?R 15.30 
'67. 164.9 00.000 OP. 070 0910 00.00 1*0000 ) .47? 

176.50 2.65 0,1.317 151.00 01.390 108.50 
034.? 06.453 28.5 07.128 ?4.5 08.899 09.155 15.85 

•68. 167.2 00.000 OR. 360 0923 00.00 1.0000 1.498 

1.84. 06 2.64 01.268 1 57.00 01 .33, 1 13.00 
29.8 06.834 ?5.7 08.508 08.773 16.30 
08.650 09?7 00.00 1.0000 1.524 

192.50 2.64 01.22? 165.00 01.288 118.00 

037.7 05.931 31.2 06.559 27.0 OR. 155 08.43? 16. R5 
'70. 173.0 00.000 08.940 0950 00.00 1.0000 1.550 

199.00 2.6? 01. IP? 170.00 01.243 122.00 
039.6 05 .686 32.7 06.331 ?8.3 07.836 08.086 17.45 

•71. 175.0 00.000 09.250 0964 00.00 1.0000 1.576 

206.00 2.63 1.1400 176.00 01.199 126.00 
041.3 05.475 34.? 06.11? 29.6 07.545 07.777 18.00 

•72. 178.6 00.000 09.560 0997 00.00 l. 0(100 1.602 

214.00 2.6? 01.100 183,00 01.155 131.00 

043.2 05.21) 35. R 05.R61 30.9 07.190 07.427 18.55 
'73. 1,80.9 00.000 09.P80 0991 00.00 1.0000 1.628 

222.00 2.6? 01,06] 190.00 01.11,4 136.00 

045.? 05.015 ?7.5 05.6?7 32,4 Oft. 898 'o7. 1 ?R 19.05 

'74. 1R3.9 00.000 10.200 1004 00.00 1.0000 1.654 

231.00 2.6) 01.025 197.00 01.075 142.00 

047.2 04.823 39.1 05.451 33.8 06.640 06*871 1,9.60 
'75. 1 86. ? 00.000 10.500 1017 00.00 1.0000 1,680 

239.00 2.61 00.990 204.00 01.037 147.00 

049.3 04.625 40.9 05.253 35.3 06.357 06.594 20.10 
'76, 190.2 00.000 10.900 1031 00.00 1.0000 1.706 

247.00 2.60 00.956 211.00 01.001 152.00 
051.5 04.443 42.6 05.060 36.9 06.106 0 6. 3 57 20.45 

•77. 19 ?.) 00.000 11.200 1044 00.00 1.0000 1.73? 

256.00 ?.60 00.923 219.00 00.067 
53.7 04.273 44.5 04.87? 38.4 05. PRO 

00.000 11.600 1057 


>.0 


00.00 1.. 

00.893 225.00 


1 57.00 

06.094 21.35 
0000 1.758 
00.934 1,61.00 
40.1 05.665 05.880 21.95 
1071 00.00 1.0000 1.784 
00.063 232.00 00.903 167.00 
8.3 03,939 4P. 3- 04,52? 41.8 05.415 05.629 22.50 
. 200.6 00.000 12.300 1084 00.00 1.0000 1.810 

231.00 2.58 00.835 240.00 00.87? 172.00 
0.8 03.779 50.4 04.349 43.6 05.280 05.413 23.10 
. 204.4 00.000 12.600 1097 00.00 1.0000 1.836 

2R9.00 2. 58 00.808 247.00 00.843 177.00 
3.5 03.632 52.6 04.201 45.4 04.998 05.200 23.80 


44.5 

195.1 

263.00 2.59 
04.114 46.4 04.710 

1,97.2 00.000 11.900 

272.00 2.59 
8.3 03,939 48.3- 


02.199 

01.845 031.01 07.178 
.18? 

02.1?) 

01.775 032.63 06.856 
.188 

02.046 

01.710 034.31 06.566 
.194 

0*1.975 

01 .649 036.06 06.29? 

.200 

01.909 

01.579 037.86 06.088 
.206 

01.845 

01 .535 039.76 05..R20 

.21 2 

01 .7B4 

01.48? 041.69 08.581 
.218 

01.726 

01.453 045.69 05.566 
.224 

01 .67? 

01 .586 045.79 OS. 1 61 
.230 

01.619 

01.541 047.9? 04.97? 
.245 

01 .570 

01 ,297 050. 1 ? 04.7S5 
.260 

01.52? 

01 .?5S 052.44 04.569 
.275 

01.476 

01.216 05^.79 04.405 
.290 

01 .177 057’. 21 04.251 
.305 

01.591 

01.140 059.70 04.064 
.320 

01. 

01 .106 062.27 05.91,5 
.335 

01.31? 

01.07? 064.9? 05.774 
.350 

01.276 

01.040 067.64 03.620 
.365 

01.241 

01.008 070.55 03.482 
.380 

0,1.207 

00.979 073.60 03.349 
.385 
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TABLE I- (Continued) 


FL ND. WGT • - VC IK)- VC <L ) J RK— -flXY— DMGAK KALPHA KRFTA- I.AI.PHA — CK 

-MK- KFDftfi CKL- NKL LI— -CL) — — L?— -CL? L3 CLM— — Ml — 

-CM l- -M? CM?- — M3 — CM3- — M4 — — M5 — C,MN — »-Nl — DMGAI. 

• pr * r? • ?n 7 .? on.ono i 3 .ono nil no.no l.oooo i.ra? , ni.175 

?9fl.nn ?. 37 00,78? 355.00 00.815 1R3.00 00.950 076.84 03.318 
086.1 03.484 5^.7 04.04? 47.3 04.797 04.994 ?A.40 .390 

' R I * 83. ?09 • 0 00.000 13.400 11?A 00.00 0.6667 l.RRR 01.J44 

307.00 ? , 56 00.757 ?6?.00 00.789 18R.00 00.9?3 079.94 O^.nqq 
068.9 03.355 57.1 03.90? 49.? 04.61? 04.808 ?5.10 .395 

*TH*90. ?3? , 1 00.000 16.300 1?17 00.00 1.0000 2.070 00.958 

00.761 109.40 03.401 

094.3 0?.577 78.1 03.080 67.5 03.573 03.745 31.30 9.4Q? .4?n 

• 11*9?. ? 3 8 . ) 00.000 17.300 1?44 00.00 1.0000 ?.!?? 00.91). 

00.7?? 1 ? 0 .70 0? . ?3 5 

10^.0 0?.394 86.? 0? . 88? 74.5 03.300 03.496 33.60 8,814 .440 


The order in which the Element File is used as input is most important. If 
we run in MODE = 1, the order of the element cards as submitted to the computer, 
must be the same as the order in which the probe measured the elements; for 
example, if data set 1 consisted of elements Ei, E 2 , E 3 , measured on spectrom- 
eters 1, 2, and 3, respectively, and set 2 consisted of E 4 , E g , measured on spec- 
trometers 1 and 2 respectively, then the order of the element cards must be E lt 
E 2 , E 3 , E 4 , E s . If other elements are present in the unknown and their weight 
percentages are known, their element cards come next. If oxygen is present in 
the unknown, i(p element cards are next. The last elements considered in the 
Element File are the file cards of elements present in the standards but not in 
the unknown. The element file is input following the DATA File. Figure 9 shows 
data for Input Files A and B for one of the sample problems (Appendix VI). ~ 

C- Composition Files 

The composition file consists of the following types of data: 

(1) Title card for the unknown sample. On this card may be punched up to 
80 characters of alphanumeric information to aid the experimenter in identifying 
his problem. This information is stored in the UTITLE array. The statement 
in MAIN; 

Read (IN, 40) (UTITLE(I), I = 1,20) accomplishes the data transmittal. 

(2) Composition card(s) for the unknown sample. The composition of each 
element in the unknown (wt. %), if known before the analysis, is input here. If it 
is not known it is set to zero. For example in columns 1-10 is punched the 
weight % of E(l) present in the unknown; in columns 11-20, the weight % of 
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E(2), . . in columns 71-80, the weight % of E(8). On the second card in columns 
1-10 goes the weight % of E(9), and so forth, until the weight percentage of each 
of the NEU elements in the unknown, if available before electron probe analysis, 
has been entered (1 ^ NEU ^ 10). Note: The order of the elements in the unknown 
should be the same as in the E list. If none of the sample compositions is known, 
they are all set to zero. The statement in MAIN 

READ (IN, 2) (C(I), 1 = 1, NEU) 

2 FORMAT (8F10.2), accomplishes the data transmittal. 

(3) Cards for the Standards. These cards are input using the following 
statements : 

DO 600 K = 1, NM 

READ (IN, 40) (STITLE(I,K), I = 1,20) 

READ (IN, 41) NEST, (ESTD(I,K), I = 1, NEST) 

READ (IN, 2) (CS(I,K), I = 1, NEST) 

600 CONTINUE 

The first READ inputs the title card for standard #K. The second READ inputs 
the second card which contains the number of elements in standard #K and the 
chemical symbols of each of these elements. On this card columns 1-4 contain 
the number of elements, 5-8, the name of ESTD (1,K), 9-12, the name of ESTD 
(2,K), . . . ., 41-44, the name of ESTD (10, K). All entries on this card must be 
right-adjusted. The third READ statement inputs the next cards which give the 

wt. % of ESTD (1, K), ESTD(NEST,K), respectively. The same format 

(8F10.2), as that used to input any known compositions in the unknown sample, 
is also used here. The order of the standards, K ® 1, 2, . . ., NM, must be the 
order in which their intensities were measured* Thus, if STANDARD #12 was 
measured first, then information reg§i#Jif B will be found whenever K = 1; if 
STANDARD #10 read second, then regarding it stored for K = 2. 

(4) Cards for the backgrounds. II K8QB > 0, then for each background 
standard a title card of identifying information is input. This is accomplished 
by the statements in MAIN: 


DO 800 K = 1, NBGS 

800 READ (IN, 40) (BTITLE(I,K), K = 1,20). 
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Figure 10 shows data from the Composition File "as used in the sample prob- 
lem illustrated in Appendix VI. 

D-Intensity Data File 

Whenever the program is operating in MODE 1, subroutine DATPTS inputs 
D. The following procedure may have to be changed according to the type of out- 
put data obtained from each particular electron probe. The type of identification 
scheme used at GSFC is outlined here. To change the scheme, DATPTS, Appen- 
dix in, may have to be reprogrammed. 

Each probe measurement is recorded on a punched card, which is divided 
into 10 fields. These fields contain 

(1) the x-ray counting time (TIME) 

(2-4) the intensity counts for each spectrometer (XINT(I)), I = 1,2,3) 

(5) the sample -type code (SAME) for unknown, standard, or background, 

(6) the set-number (NUMRDG), which tells whether element group 
1 (E(l) - E(3)), 2 (E(4)-E(6)) or 3 (E(7)-E(9)) is being measured 

(7) the group - number (INDIC) which tells the number of times the standards, 
as a group, have been measured. INDIC is used so that the drift calcu- 
lation can be made several times during a run 

(8) the number or identification of the standard or background being measured 
(NOST). As discussed in the description of the Composition File, the 
background and standards are remembered according to the order in 
which their intensities were measured (1,2,3. . . .). Therefore this re- 
alignment is used for NOST 

(9) the dummy variable (IUNK) for unspecified usage, 

(10) the end of data card (LAST). This field is used on the last data card of 
a problem to indicate the end of the data. Its usage allows a second 
problem to be run. 

The statements: 

READ(IN,1) TIME, (XINT(J), J = 1,3), SAME, NUMRIDG, INDIC, NOST, IUNK, LAST 
1 FORMAT (F6.0, 3F8.0, 411, 12,11) 

input one card. The data are then processed according to the value of SAME. 

Table II presents a capsule review of the use of the 10 fields with some additional 
comments. 
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TABLE II 

Format for Intensity Data File 


Field 

Columns 

Format 

Parameter 

Comment 

1 

1-6 

F6.0 


duration of measurement 

2 

7-14 

F8.0 


intensity in counts from spectrometer #1 

3 

15-22 

F8.0 

n 

intensity from spectrometer #2 

4 

23-30 

F8.0 

XINT(3) 

intensity from spectrometer #3 

5 

31 

11 

SAME 

SAME = 1 unknown measured 
SAME = 2 standard measured 
SAME = 3 background measured 

6 

32 

11 

NUMRDG 

NUMRDG - 1 measurements taken on 
E(l),E(2),Ec3) (Set 1) 

NUMRDG = 2 measurements taken on 
E(4),E(5),E(6) (Set 2) 

NUMRDG = 3 measurements taken on 
E(7),E(8),E(9) (Set 3) 

7 

33 

11 

INDIC 

INDIC = i, i ~ 1,2,. . . ., 9 

ith time readings are taken on standards 

8 

34 

11 

NOST 

NOST = i, i = 1,2,. . NM measurement 
taken on the ith standard 

9 

35-36 

12 

IUNK 

IUNK = i, i = 1, 2, . , , ., 99 
Dummy variable for other usage 

10 

37 

11 

LAST 

LAST = 0 except on last data card where 
LAST = 1. This indicates the end of 
the data for the problem. 

i 
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APPENDIX V 
Output from the Program 


The format of the cL\ta output and the FORTRAN statements which control 
this output are given below. The FORTRAN statements are referred to by their 
sequence numbers (SN) which appear to the right of each statement in the pro- 
gram listing (Appendix I). 

Data set B, that is the E file, and selected parameters from data set A are 
printed out by MAIN, SN 15200: 19700. The matrix of mass absorption coeffi- 
cients is calculated in SUBROUTINE ABSCO, SN 6500: 8100. The raw probe 
data are output in SUBROUTINE DATPTS at SN 4700: 4900, SN 6200: 6300, SN 
7100: 7200, according to whether the measurement was of an unknown, a standard, 
or a background, respectively. For each standard, the title and composition 
cards are printed by MAIN, SN 24100: 24300; and the atomic number correction 
is also printed by MAIN, SN 24400: 24400. For each element in each standard, 
the absorption correction is printed by MAIN, SN 25700: 25700 and the fluorescence 
correction is output in FLUCOR, SN 5200: 5300, SN 5800: 6200. 

Whenever a new standard combination ;is used, a message is written by 
MAIN, SN 27500: 27500 telling which standard is used for each measured element 
in the unknown. The sample compositions are printed at each iteration step: 

MAIN, SN 32300: 32600 only for the first set of intensity data. The weight % sum 
is also printed out. No attempt is made to normalize the sum to 100%. 

After the final iteration on the first data point in a MODE 1 problem and for 
MODE 2 problems, the atomic number corrections, the absorption corrections 
and the fluorescence corrections are printed by statements MAIN, SN 31000: 

31600, FLUCOR, SN 5200: 5300, FLUCOR, SN 5800: 6200 respectively. For sub- 
sequent data points the compositions of the first and last iterations are written, 
MAIN 29200 and 32600. For MODE 2 problems intensity ratios are printed by 
MAIN, SN 33300: 33300. Sample output data are given in Appendix VI, where 
several sample problems are considered. The output format can be easily 
changed by altering the FORTRAN statements previously indicated. 
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APPENDIX VI 
Sample Problems 

Three sample problems will be described in order to illustrate the use of 
this program. A copy of the input data (Files A-D) and the output for each of 
the problems is supplied as a guide for the successful application of the 
program. 

Problem 1 - (Geological Application)* A sample is analyzed in MODE 1 
for the following 8 elements: Ca, Mg, Si, Al, Na, Mn, Cr, and Fe. Oxygen is also 
present as is 0.1 wt % Ti. Six standards are used: pure diopside for Ca, Mg 
and Si, and 5 other standards, one for Al, Na, Mn, Cr and Fe. The data were 
taken on an electron probe with 3 spectrometers and a takeoff angle of 52.5°. 

The operating voltage was 20 KV; only K lines were measured. Input Files 
A-C have already been illustrated (Figures 8-10). Input File D, the Intensity 
Data File, is shown in Figure 11. A shortened version of the Output Data is 
given in Figure 12. 

The initial estimate of the unknown composition using a linear relation 
(Equation 15) is given in Table III as well as the final calculated compositions 
determined after 3 iterations. 


TABLE III 

Calculated Compositions Pm.-* \l ± 


Element 

Initial Estimate (wt. %) 

Final calculated composition (wt. %) 

Ca 

11.54 

11.57 

Mg 

12.55 

12.87 

Si 

25.13 

26.1 

Al 

1.15 

1.08 

Na 

0.88 

0.93 

Mn 

0.09 

0.09 

Cr 

0.59 

0.58 

Fe 

2.57 

2.59 

Ti (known) 

0.10 

0.10 

0 

43.78 

45.0 

Sum 

98.4 

100.9 


For all elements measured in this problem, the only correction factor of 
any significance was the absorption correction. 


*This problem was kindly loaded to us by Dr. J. Boyd of the Geophysical Lab,, Washington, D.C. 


48 










Problem 2 - (Metallurgical Application)* A binary sample is analyzed in 
MODE 1 for Ti and Nb only. Pure standards of the elements are used. The data 
was taken on an electron probe with 2 spectrometers and a takeoff angle of 52.5°, 
The operating voltage was 20 KV and Ti-K and Nb-L radiations were measured. 
Input Files A-D are illustrated in Figure 13 and the Output Data is given in 
Figure 14. 

The initial estimate of the unknown composition using a linear relation 
(Equation 15) is given in Table IV as well as the final calculated composition 
determined after 3 iterations. 


TABLE IV 

Calculated Composition Problem 2 


Element 

Initial Estimates (wt. %) 

Final Calculated Composition (wt. %) 

Ti 

32.7 

34.8 

Nb 

62.2 

65.5 

Sum 

94.9 

100.3 


For Ti, an atomic number correction of +7.6% is necessary while for Nb a 
negative atomic number correction of 5% is calculated. A significant absorption 
correction (~13%) was calculated for the Ti. The absorption correction for Nb 
was insignificant. 

Problem 3 — (Metallurgical Application) A sample containing 4 elements has 
the following composition Al-lO wt. %, Mg-50 wt. %, Au-40 wt. %, Cu-0.10 wt. %. 

It is of interest to calculate using MODE 2 the expected concentration ratios 
(I /1(A)) for the 4 elements if pure standards of each element are used. The 
data is calculated for an electron probe with a takeoff angle of 52.5° and an op- 
erating voltage of 20 KV. The following x-ray linus were considered Al-K, Mg-K, 
Au-L and Cu-K. Input Files A-C are illustrated in Figure 15 and a shortened 
version of the Output Data is given in Figure 16. 

The calculated intensity ratios are given in Table V. 

This particular combination of elements exhibits a large atomic number ef- 
fect (up to 30%) and a large absorption effect (up to 50%). Only Cu is noticeably 
fluoresced, namely by the Au-L line. Large corrections of this order, com- 
pounded with f( x ) values less than 0.8 for A1 and Mg, argue quite strongly for 
the use of intermediate, 4 element, standards. 


*This problem was kindly loaned to us by D, Beaman of the Dow Chemical Company, Midland, 
Michigan. ft 
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TABLE V 

Calculated Intensity Ratios, Problem 3 


Eloment 

Composition 
(wt. %) 

ANC * 
ATN/ATD 

ABS s 
ABN/ABD 

FLU * 
FLN/FLD 

Intensity Ratio, 
Sample to Standard = 
ANC*ABS*FLU , C(U) 

A1 

10 

1.14 

0.48 

1.0 

5.5 

X 10 “ 2 

Mg 

50 

1.17 

0.69 

1.0 

41. 

X 10~ 2 

Au 

40 

0.71 

1.01 

1.0 

29. 

X 10- 2 

CU 

0.1 

0.95 

0.98 

1.07 

.101 X 10” 2 
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FIGURE 1 : MAIN PROGRAM (CONTINUED) 
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Figure 9— Input Data — Sample Problem 1 — DATA and Element File 
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Figure 11~lnput Data - Sample Problem 1 - Intensity Data File 
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Figure 12— Output Data — Sample Problem 1 
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Figure 12— Output Data — Sample Problem 1 (Continued) 
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Figure 12-Output Data - Sample Problem 1 (Continued) 
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n'Sut Data — Sample Problem 2 
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Figure 14-Output Data - Sample Problem 2 
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Figure 14- Output Data - Sample Problem 2 (continued) 
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Figure 16— Output Data — Sample Problem 3 
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Figure 16-Output Data - Sample Problem 3 (Continued) 
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